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ABSTRACT 


This  thesis  is  an  application  Df  parameter  sensitivity 
analysis  to  aerodynamic  parameters  of  a  Bank-to-Turn 
missile.  In  the  development  a  brief  review  of  trajectory 
sensitivity  theory  is  presented.  A  linear  analysis  is 
performed  using  an  Uncoupled  Pitch  Channel  Autopilot  and  a 
Coupled  Boll-Yaw  Channel  Autopilot  of  the  missile  taken  as 
model.  Finally, a  nonlinear  analysis  is  given  for  the  system. 
Comparisons  between  the  linear  and  nonlinear  cases  are 
outlined. 
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I-    PARAMETERS    SENSITIVITY   ANALYSIS    OF    A    BANK-TO-r PEN    MISSILE 

The  determination  of  changes  ia  system  performance  due 
to  changes  in  parameters  is  of  great  importance  in  engi- 
neering analysis  and  design.  Sensitivity  questions  arise 
when  model  uncertainty  is  present ,  or  when  a  range  of  oper- 
ating conditions  is  contemplated. 
The  questions  of  parameters  sensitivity  particularly  arise 
in  the  fields  of  engineering  where  models  are  idealized, 
inexactly  identified,  or  the  systems  themselves  are  subject 
to  unpredictable  changes  with  tiae  due  to  enviromental, 
material  property  or  operational  influences  so  that  there  is 
always  a  discrepancy  between  the  physical  reality  and  the 
mathematical  model. 

Sensitivity  analysis  provides   the   engineer   with    methods    for 
investigating  or   minimizing    such    parameter   deviations. 

In  general,      the    diagram  of     a  system  can   be   represented 
by   a    single  block  as   given    on   Pig.    1.1. 


Figure   1.1        General  Representation   of   a   Dynamic  System. 

From  a  mathematical  point  of  view,  what  we  call  a 
system  is  the  explicity  or  implicity  given  relationship 
between  the  input  signal  u  and  output  signal  y.  In  general, 
u  and  y  can  be  vectors.  The  character  of  this  relationship 
is   commonly   called   the  structure   of   the   system. 
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For      example,  the        structure      of      the      system        may      be 

characterized  by  the  order  of  a  differential  or  difference 
equation,  linearity  or  nonlinearity,  the  order  of  the  numer- 
ator and  denominator  of  a  rational  transfer  function  and  the 
quantitative  properties   of    the   system   parameters. 

Typical  parameters       are      initial        conditions,  time- 

invariant  or  time-variant  coefficients,  natural  frequencies 
and   samplinq  periods. 

The  chanqe  of  the  state  or  the  chanqe  of  the  output 
variable  with  time,  can  be  caused  by:  (1)  the  influence  of 
input  signals,  (2)  the  chanqe  of  parameters.  These  quanti- 
ties   are   shown    in  Fiq.    1.2. 


2Ci) 


input 

Sipncxls 


•z     Pa-i-ameter 
^t  ChQ.nggfc 


SYSTEM 


Figure   1.2       Quantities   Affecting  the  Dynamics  of    a   System* 

This  work  addresses  the  application  of  the  sensitivity 
analysis  to  the  variation  of  aerodynamic  parameters  of  a 
Bank-to-Turn  missile. 

PARAMETER  SENSITIVITY  is  the  effect  of  parameter  changes  on 
the  dynamics  of  a  system,  say,  the  time  response,  the  state, 
or   any  other  quantity  characterizing    the   system   dynamics. 

In  order  to  have  a  realistic  model  to  work  with  a  NASA 
Contractor  Report  [ Ref •  2]  was  adopted  as  a  reference  for 
application  of  the  sensitivity  analysis  on  the  variation  of 
aerodynamic  parameters  of  a  Bank-to-Turn  missile.  Appendix 
A,  B  and  C  give  the  detailed  aerodynamic  assumption  of  the 
system. 
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One  brief  explanation  of  the  sensitivity  theory  and 
analysis  is  given  as  a  guideline  to  better  understand  the 
subsequent   work. 

An  analysis  of  the  sensitivity  of  the  aerodynamic  parame- 
ters applied  to  the  case  of  the  linear  uncouplsd  pitch  and 
coupled   roll-yaw  autopilot    is   performed. 

The  pitch  and  roll-yaw  autopilots  are  couplsd  as  given 
in  Fig.C.1  in  the  appendix  C.  A  nonlinear  parameter  sensi- 
tivity analysis    is    performed. 

Conclusions  about  the  influence  of  parameters  of  concern 
are  given  and  a  comparison  between  the  linear  and  nonlinear 
case  is  done.  Comments  and  recommendations  for  the  future 
are  outlined  in  order  to  delineate  the  continuity  of  the 
present   work. 
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II.    SENSITIVITY    THEORY 

A.  INTRODUCTION 

The  basis  of  all  sensitivity  considerations  in  the  case 
of  time-invariant  parameter  variations  is  tha  so-called 
sensitivity  functions.  Dynamic  systems  can  be  characterized 
in  several  ways: in  the  time  domain,  in  the  frequency  domain, 
or  in  terms  of  a  perfomance  index.  There  is  evidently  an 
adequate  number  of  ways  to  define  tha  sensitivity  functions 
of  a  dynamic  system.  The  definition  that  is  actually  used 
depends  on  the  form  of  the  mathematical  model  as  well  as  on 
the   purpose   under  consideration. 

The  sensitivity  functions  can  be  classified  into  the 
following  three    categories: 

(1)  Sensitivity    functions   in   the    time   domain 

(2)  Sensitivity    functions   in   the    frequency  domain,    and 

(3)  Perfomance- index   sensitivity 

In  this  chapter  we  will  outline  just  the  sensitivity  func- 
tions in  the  time-domain  of  the  continuous  systems, 
^plications  of   this    analysis  follow   throughout  this   work. 

The  reader  can  find  detailed  informations  about  the 
others  two  categories  of  sensitivity  functions  in  the 
[Ref.    1]. 

B.  TRAJECTORY  SENSITIVITY  FUNCTION  OP  CONTINUOUS  SYSTEMS 

A  continuous,  possibly  nonlinear  system  of  nth  order 
can,  in  general,  be  described  in  tha  state  space  by  a  vector 
differential  equation  as   seen   in   Eqn.2.1. 
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X  =  f(X,t,u,*o)  ,  X(to)  =  XO  (2.1) 

Here  X  is  an  n  x  1  state  vector,  f  an  n  x  1  rector  func- 
tion, u"  an  input  vector,  ©<0  a  nominal  r  x  1  parameter 
vector,  and  X  is  the  n  x  1  initial  condition  vector  or 
initial  state.  Egn.2.1  is  called  the  NOMINAL  STATE 
EQUATION. 

Assuming  that  the  parameter  vector  deviatas  from  the 
nominal  value  <*0  by  ^o^,  we  have  the  Eqn.2.2. 


X  =  f(X,t,u,«),  X(t0)  =  XO  (2.2) 


with  the  initial  conditions  X°  unchanged.  This  equation 
is    called  the   ACTUAL    STATE    EQUATION. 

Now   it    is   assumed  that    Eqn.2.2  has   a   unique   solution,    X= 

X(t,»c)      for  all  admissible    initial   conditions  and    parameter 

values. 

X  is  of  course   a  function  of  u  ,    X°  and   t  as  well. 

o 

However  ,  this  dependence  is  not  needed  for  the  following 
considerations  and  will,  therefore,  be  dropped  for  ease  of 
notation.  Furthermore,  the  solution  X  is  assumed  to  be  a 
bounded   continuous    function    in  t   and    ©<,    . 

If  the  parameter  takes  on  its  NOMINAL  valus  <x0  r  the 
nominal   solution   X0     =  Xft,©^)      is   obtained.  If,      on   the 

other  hand  the  ACTUAL  solution  is  given  by  X  =  X(t,  <K  )  , 
then  the  parameter- induced  change  of  the  stata  vector  is 
given   by 


AX(t,*)    tX(t,  X)    -    X(t,<xo)  (2.3) 
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A    first-order      approximation   of  A  X      can    be      obtained   by 
using   Taylor  expansion  in   the    form   of    Eqn.2.4. 


4I(t,«|    =     Zj    -zz—T    I-  -J  (2.») 


Ki  T>*. 


J 


This  equation  can  be  viewed  as  a  definition  of  the 
parameter-induced  trajectory   deviation. 

The  subscript  oto  shall  indicate  that  the  partial  deriv- 
ative   expressed    by  "D       is  taken   at   nominal    parameter   values. 

Let  the  state  X  of  a  continuous  system  be  a  continuous 
function  of  a  time-invariant  parameter  vector  o<  - 
{ <x,  f^,...^,.}.  Then  the  partial  derivative  can  be 
defined   as 


■& 


""DUjSt) 


Vet,**)  ■  -a-v  «*  <2-5> 


W#    2,    .    .   .   ,    n 

j— 1  ,     2,0     •     •    •     i     T 
fiqn.2.5  is      called     the  trajectory     sensitivity     vector      with 
respect   to   the    jth    parameter. 

Note  that  the  trajectory  sensitivity  vector  is  of  the 
same  dimension  as  the  state  vector,  namely,  n.  Its  compo- 
nents   are   the    TRAJECTORY    SENSITIVITY    FUNCTIONS    as 


Eqn.2.6   is    the    partial    derivative   of   the   ith   state   vari- 
able   in  relation  to    the   jth    parameter   as 
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(2.7) 


Hence,  all  n  x  r  trajectory  sensitivity  functions  form 
the  trajectory  sensitivity  matrix  as  given  in  Eqn.2.8  or 
2.9. 


A       -    {  A.  ,..  •  ,  >,„.}     - 


(2.8) 


"  D  X 


A: 


1)* 


*  »  • 


~a  xi 


T)Xn  lxn 


Hot 


^>^»-   . 


(2.9) 


The  columns  of  A  are  the  trajectory  sensitivity 
vectors  X^  .  Here  X  is  the  Jacobian  matrix  3f  the  state 
vector  with  respect  to  the  parameter  vector  oc  ,  taken  at 
the    nominal   parameter  values. 

With  these  definitions  the  parameter-induced  change  of 
the   trajectory   can   be  taken    as 


A  x(t,°c)   =  A   (t,^)**  »  ^P    X;  A<*~j 


(2. 10) 


J-' 


Where  o<    =  <*0  ♦   A  *<  ,        which      is    the      ACTUAL      parameter 
vector  of  the  system. 
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C.       TRAJECTORY    SENSITIVITY    EQUATIONS    OF    CONTINUOUS    SYSTEMS 

Lets   consider  the  general   continuous  system   described  as 
previously   by  the   state   equation    (Eqn.2.1). 


X    =    f  (X,tfu,<*),    X(t0)     =    XO  (2.11) 


Where  X  denotes  the  n- dimensional  state  vector,  c*  the 
r-dimensional  parameter  vector,  f  an  n-dimensional  vector 
function,  and  u  the  input  vector  independent  of  c<  .  it  is 
assumed  that  the  continuity  conditions  are  fulfilled  and 
that     «.     is  time-invariant. 

Considering  «  -parameters  and  taking  the  partial 
derivative  of  X  (Egn.2.11)  with  respect  to  °<j  ,  one  obtains, 
by   the  application   of  the  chain  ruls 


^X  j)4       ^x  ~M  ^X° 

^•y      ~^*     ^**o      ^*j       ^W 


(2.12) 


The  derivative  of  the  initial  conditions  vector  X°  with 
respect   to    OCj    is  zero,    since   X<>    doss   not   depend    on     o<    . 

If  X  is  r-dimensional,  there  are  r  equations  of  the 
form  of  Eqn.2.1 2.  If  we  now  interchange  the  seguence  of 
taking  the  derivative  with  respect  to  time  t  and  **j*  ,  and 
then  let  <x.  approach  to  °*o  t  ona  obtains  ths  following 
equation 


Aj    = ™L\'4    ZZ L'    A^)-    0  (2.13) 
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Here  A;  =<-^. —  -r-  |ZT  is  the  trajectory  sensitivity  vector 
with    respect  to    the    jth    parameter. 

The  Eqn.2. 13  is  called  the  sensitivity  equation  in  the 
state   space   or    the    TRAJECTORY    SENSITIVITY   EQUATION. 

The  abovs  shows  that  for  c<  -parameters  all  initial 
conditions  of  the  trajectory  sensitivity  equations  are  equal 
to   zero. 

Now   consider  the   output    vector  aquation,    as  given   by 

y      =7(Xrt,u,X)  (2.14) 

In  a  procedure  similar  to  that  above  one  obtains  the 
alqebraic  sensitivity  equation  as  seen  in  Eqn.2. 15. 

Cj-=-- ,  «  ------ 1_  Xj  ♦  ------ 1--  (2-15) 


<* 


_    _  Which     relates        the     output        sensitivity        vector  yj"   ■ 

o^  T  . 

r-c-T  —  |—    to  the   trajectory      sensitivity   vector    Aj     .        This 

equation   is   called   the   VECTOR   OUTPUT    SENSITIVITY    EQUATION. 

Using  the  trajectory  sensitivity  matrix  A  and  the 
output  sensitivity  matrix  CT*  ,  the  above  result  can  be 
rewritten  in  the   following    general  form 

X      *   -— — U  A    ♦    — j  — 1_       r\o   =    o  (2.16) 


CT      =    ----- UA    ♦    —£--,_  (2.17) 

6>  X       *o  (Joi        °<o 


28 


•These  equations  are  called  the  STATE  SENSITIVITY 
EQUATIONS  of  the  system.  It  is  seen  that  these  equations 
are  linear  whether  the  original  system  is  linear  or  nonlin- 
ear. 

If,  in  particular,  the  original  system  is  linear,  the 
state   equations,    take  the  form 


X    -    A    X   «■    B    u     ,    X(to)     =   XO  (2.18) 


Y    =    C    X   ♦   D    u  (2.  19) 


Bhere,    in  general, 

A~=I(o<)     ,  "B   =    ?(<*)     ,    C    =    C  (<*)     ,D=r7(*<.), 
~X  =    X~(t,o<.  )    ,    amd  7  =  7(t,&<  )  . 

Note,  however,  that  u  is  not  a  function  of  <k  if  u  is 
defined  as   an  external  input  of  the   system. 

Now  taking  the  partial  derivatives  with  respsct  to  °<J  , 
reversing  the  order  of  differentiations  with  respect  to 
time  and  c^ij  ,  and  letting  ££  approach  to  ©<a  ,  the 
TRAJECTORY    SENSITIVITY  EQUATIONS    are    obtained 


X,    *   A      V  *    -  - 1    XO   ♦    --— -|_    U(t)     ,  Att0)  =  G    (2.20) 


Where  A0  =    A(<*0)    ,    x^  ■  X(t,  <*Q)  ,   and   j=    1,      2,...,     r. 

initial     condition  vector 
X(t0)    does   not    depend  on   ©<    . 


The    initial     condition  vector    X -Lto)*      is   again      zero   since 
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By  similar      procedure  the     VECTOR   SENSITIVITY        EQUATION 
becomes 


(2.21) 


Where  C0  =    C ( c^ )    and  j=1,2,...,    r. 

In  the  case  of  linear  systems,  the  vector  sensitivity 
equations  have  the  same  A  matrix  as  the  nominal  state  equa- 
tions and   hence   the    same   characteristic    matrix  si    -    A. 

They  differ  from  the  nominal  original  state  equation  only 
in  the  driving  function  and  the  initial  conditions.  The 
latter  are  all  zero.  The  driving  functions  can  be  obtained 
by   solving   the    nominal  state  equations. 

Here,  in  the  linear  case,  the  sensitivity  e?uations  can 
be  joined  to  the  system  equations  ,  forming  the  so-called 
COMBINED   SYSTEM    as    given   in    Eqns.2.22    and    2.23. 
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The  simultaneous     solution   of     these   differential      equa- 
tions    using  the     same  standard      procedure   for     each   of      the 
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above  matrix  equations  yields  the  output  vector,  the  state 
vector,  and  the  corresponding  output  and  trajectory  sensi- 
tivity vectors  Gy  ,  and  Aj  •  If  there  are  r  parameter 
variations,  r  systems  of  equations  of  the  above  type  have  to 
be  solved.  Since  the  homogeneous  part  of  the  original 
differential  equation  is  identical  with  the  homogeneous  part 
of   the  sensitivity    equations   with   rsspect   to   all    parameters. 

A  graphical  interpretation  of  Eqns.2.22  and  2.23  is 
given    in  Fig. 2- 1 . 

Since  the  driving  function  of  the  sensitivity  model 
contains  the  nominal  state,  the  measuring  circuit  for  the 
trajectory  sensitivity  functions  consists  of  a  connection  of 
the  nominal  original  system  and  the  sensitivity  model  as 
illustrated  by   Fig. 2.1. 

If  the  actual  input  u  is  applied  to  this  structure  the 
trajectory  and  output  sensitivity  vectors  Xj  and  CV  can  be 
measured  at  the    points   1   and  2   respsctively   in  the   Fig. 2.1. 

In  order  to  measure  all  r  sensitivity  vectors  simultane- 
ously,  r  sensitivity   models    are  required. 

D.       STRUCTURAL    METHOD 

The  basis  for  the  measurement  of  the  trajectory  sensi- 
tivity functions  is  the  STRUCTURAL  INTERPRETATION  of  the 
trajectory   sensitivity  equation. 

The  physical  system  represented  by  the  trajectory 
sensitivity  equations  is  called  the  trajectory  sensitivity 
model   or  the  sensitivity  model   in   the   stare   spaca. 

Regardless  of  the  nature  of  the  original  system  the 
sensitivity  model  in  the  state  space  is  always  linear.  The 
graphical  illustration  of  sensitivity  equation  can  be  given 
by  Fig. 2. 2.  The  system  described  by  the  sensitivity  equation 
is  refered  to  as  the  sensitivity  model  of  the  original 
system.      For  each  output,      a      system    has   as   many    sensitivity 
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models  as  parameters  of  interest.  Both  the  nominal  original 
system  and  the  corresponding  sensitivity  models  form  the 
COMBINED  SYSTEM.    The   sensitivity   model   is   always    linear.      If 

Y(t,<*.)    is   the    output  of   a    system,      the    corresponding   sensi- 

—  —  (J  V 

tivity   functions      5Ty  (t,  oCp  )      =  "p<Vo4:~t   are  the      outputs   of 

the  corresponding  sensitivity  models.  Consequently,  in  order 
to  measure  the  output  sensitivity  functions  simultane- 
ously,the  nominal  original  system  and  the  sensitivity  models 
have  to  be  measured  at  the  outputs  of  the  sensitivity 
models. 

The  double  frame  used  for  the  original  system  in  Fig.  2.2 
is  to  indicate  that,  in  general,  the  system  may  be  nonlin- 
ear, whereas  the  sensitivity  models  are  always  linear.  In 
the  nonlinear  systems  the  Egn. 2.12  applies  and  different 
programs  have  to  be  set  up  for  the  original  and  the  sensi- 
tivity equations. 

For  some  applications,  the  sensitivity  functions  with 
respect  to  all  parameters  are  required  simultaneously.  If 
there  are  r  parameters  of  interest,  r  sensitivity  models 
would  be  needed.  This  implies  a  rather  extensive  computer 
time. 

A  method  of  determining  all  output  sensitivity  functions 
of  a  system  simultaneously  by  a  single  sensitivity  model  is 
available  ,  which  is  also  called  the  method  of  SENSITIVITY 
POINTS.  This  method  has  application  just  for  linear  systems. 
Detailed  explanation  about  the  application  of  the  sensi- 
tivity points   theory   can  be    found  in    the  [Ref.    1]. 

E.       OUTPUT    SENSITIVITY   FUNCTION    OF   CONTINUOUS    SYSTEMS 

Consider  the  input-output  behavior  of  a  continuous, 
possibly  nonlinear,  single-variabla  system  described  by  an 
ordinary   differential  equation   of   the   type 
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f{   ym)  ,y(n",;    ,...,y,t,   <x0  }    =   0  (2.24) 


With  the  initial  condition  y(t0)  =  y.  °,  i=0,  1,...n-1.  y 
denotes  the  output  signal,  t  the  time,  and  <x  a  single 
time-invariant      or      slowly      varying  -parameter      that      has 

nominal   value    «^o    . 

In  general,  f  is  a  function  of  the  input  u  as  well. 
However,  if  u  is  an  external  input  which  does  not  depend 
on  <x„  ,  the  dependence  of  f  on  u  is  not  relevant  in  further 
considerations. 

Let  us  suppose  that  ths  above  NOMINAL  differential  equa- 
tion   has  the  unique    solution  given  below 


*o    s   Y(t,  *o>  <2-25> 


Which  one  shall   call   the   NOMINAL    SOLUTION. 

Let  us  now  assume  that  the  parameter  changes 
from  X0  to  <X>  =  o<0  ♦  A^,  where  ^  ©<w  is  time-invariant  or 
slowly  varying  with  time.  <K  is  called  the  actual  parameter 
value.  The  corresponding  ACTUAL  DIFFERENTIAL  EQUATION  can, 
then,    be   written   as 


m)     m-»; 
f    =    {3    ,y  ,...,y,t,«}    =   0  (2.26) 


Note     that    by     this     change  of    <Xq    into    oc      the      initial 
conditions   remain  unchanged,    namely   y  (t. )    ■   y.°. 
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The    corresponding  solution    is   given   as 


y  =  y(t,<x)  (2.27) 


Which   we  shall    call    the    ACTUAL    (or  PERTURBED)    SOLUTION. 

It  is  assumed  that  y(t,<x)  is  of  the  same  type  as 
y(tf  «<„),  and  y(t,  x)  deviates  infinitesimally  from  y(t,°<o) 
if  <X  deviates  infinitesimally  from  °<o  •  The  conditions 
for  fulfilling  this  requirement  are  given  in  the  mathemat- 
ical literature.  For  our  purpose,  it  is  sufficient  to  know 
that  y  is  continuous  in  ^  if  f  is  continuous  in  y  which  is 
true    for   all  continuous  systems  and     °<o    +   0- 

With  the  above  assumptions,  the  actual  solution  y  (t,  <X0 
♦  Zot  )  can  be  expanded  into  a  Taylor  series  around  °<0  , 
yielding  the  Egn.    2.28. 


y(t,*)   =  y(t,*g   ♦^™L*ot  +  -5"~.T£'  ***  *•••  (2.28) 


If  A<x  <<  X0  ,      the  Taylor   series    can   be  truncated  at  the 
linear  term. 
This   gives  the    Egn. 2.29. 


y(t,<x)    =    y(t,0<o)    ♦    -- |     ^<*  (2.29) 

7)  <*     <*o 


For     finite      values     of  A  c<  ,      this     expression     can     be 
considered   a  first-order  approximation   of   y(t,<X.). 
The   actual   output  can  be  written    as 
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y(t,cx)    =    y(t,*0)     ♦    <Mt,<XD)    adc  (2.30) 

Where  y(t,^0)  is  the  nominal  output  and  <T  (t,  <*0)  is 
the    parameter-induced  output  error. 

With  0  (t,  ©c0)  defined  as  -"^"Jj — H^  •  The  parameter 
induced   output    error   is   in    these    tarms   defined  as 


(t,*  )    =     6*  (t,  «o)    Aoc  (2.31) 
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Figure   2.1        Graphical  Interpretation    of    Eqns.2.2  2   and   2.23. 
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Figure  2.2    Strutural  Diagram  of  the  Combined  System. 
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III.    APPLICATION   OF    SENSITIVITY    AN£IiXS.IS_   TO    LINEAR    SYSTEMS 

A.  INTRODOCTION 

The  sensitivity  theory  described  previously  will  be 
applied  to  the  case  of  an  uncoupled  pitch  autopilot  and  to 
a  roll-yaw  coupled  autopilot  presented  in  the  appendix  C. 
For  each  case,  the  linear  equations  of  the  nominal  system 
are  presented  in  state  variable  form  and  the  correspondent 
trajectory  sensitivity  equations  are.  derived.  A  sensitivity 
analysis  of  the    systems   is    performed. 

B.  UNCOUPLED  PITCH  AUTOPILOT  ANALYSIS 

1.      Linear   Equations  of    the   Nominal   System 

From  the  block  diagram  of  the  Fig.B.1  in  the 
appendix  B  one  can  obtain  the  following  nominal  state  equa- 
tions of  the  uncoupled  pitch  autopilot. 


X1    =  C2    A3  X2    +   C2    A4    X3  (3.1) 


X2    a   X1    -    K    C1    A1    X2  -    K    C1    A2    X3  (3.2) 


X3    =   -  C3   X3    ♦  C3    Conv   X6  (3.3) 


X4    =    -   CI   CH    A1    X2    -   C1    C4    A2    X3   -C4    X4  (3.4) 
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X5  =  C7  Xft  -  C5  X5  -  C6  NZC 


(3.5) 


X6  ■  C9/Conv  ♦  C2  C8/conv  A3  X2 


(3.6) 


♦  C2    C8/Conv    Aft    X3   -   C7    C8    Xft    ♦    (C5    C8    -    C9)     X5 

♦  C6    C8    NZC 

For  the      purpose   of      using  the      linear   method      these 
state   equations    must   be   presented   in    matrix   form   as   given   by 


0  C2A3  Q2*A 

1  -K1CfAt      -KCtA2 


O  o  -c3 


o     -c1qqA1     -ciS*z 
o  o  o 

c 


o 
o 
o 


o 
o 
o 
-c- 


,-c 


conv       qonv        co«v 


0 
O 
CjConv 
O 
O 
O 


jx,: 

- 

o- 

y2 

o 

*j 

o 

*4 

+ 

o 

*s 

-% 

*6 

M 

N; 


(3.7) 


-«    u    — 


The  correspondence    of   the   state   vectors   is: 

X1    =  Cj     ,    12    s  «(    ,    13    a    Jp  ,     X<*    =    X    ,     X5    =    Y    ,    X6    =    Sp 

Definition  of  the  constants  C1  through  3  9  are  given 
in  appendix  B.  The  parameters  of  interest  for  this  system 
are   given  by  11,   A2r    A3,    and  Aft   which   are   : 

A1    =    C..      ,    A2    =    CM^        ,    A3    =   Cm       ,    Aft    =    C         . 

This  nomenclature  is  used  here  to  easily  apply  the 
sensitivity  theory  and  to  avoid  the  inconvenience  of 
carrying  symbols  and   constants. 
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2.      Sensitivity    Equations 

To      apply     the      procedure   developed      in      chapter     2, 
consider   the  trajectory  sensitivity   equation    (Eqa.2.13). 


X;    =    kQ    X:    ♦     — l_XQ    ♦ 1     U(t),    A;^}:0       (3.8) 


From  Eqn.3.7   one   can   see   that 


V 


o 

1 

o 


C2A5 


C2A4 


-K1CfAt      -KCfA2 


-C, 


O 
O 
O 


o     -c^a,     -Sq,^     -c4 


o 
o 
o 

-C- 


.0 
o 
Cjconv 
o 
o 


Cg_    Cgce  A3    C&C6^      ^  C-C.-C0 


cony       conv 

u     — 


comv 


(3.9) 


B:[o       O 


-  c 


c 


c6  Ca  3 


(3.10) 


The  partial  derivatives  "■^w""' '  ^  an  d  ~J"(SSf~'o4~are 
evaluated  considering  the  parameters  of  interest  A1,  A2r  A3 
and  A4,  that  are  respectively  the  aerodynamic  coefficients 
present   in   the    pitch   channel. 
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Applying   the    partial  derivatives   with   respect   to  the 
parameters    we  have   the  following   matrices 


DA. 


O 
O 

o 
o 
o 
o 


o 
o 
o 
o 
o 
o 


C2  o 

o  o 

o         o 


o 

CiC9 
Cor\V 


o 
o 
o 


o 
o 

6 
o 
o 
o 


o 
o 
o 
o 
o 
o 


(3.11) 


TA_ 
"&a" 


o  c2       o  o 

o  o        o  o 

o  o        o  o 

O  O           O  O 

o  o        o  o 

o  £^2    o  o 
Conv 


o  o 

o  o 

o  a 

o  o 

o  o 

o  o 


(3.12) 


"DA. 

DA3 


o 

o 
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o      o 


o 
o 
o 
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o  o 

o  o 

o  o 

o  o 

o  o 


(3.13) 


"&a. 


O     -KCj     O 


o  -c, 
o  o 
o      o 


,Qq  o 

o 
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o 
o 
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o 
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o 
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o  o 

o  o 

o  o 

o  o  1 


(3.14) 
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and 

because   B   is  independent   of    the   parameters   of   interest. 
In   terms  of    components  one   can   see 


L  J  "ft  tfj  %j 


(3.15) 


i=1,2,...,6    j-1,2,3    and  4 

Here,  four  models  are  necessary  to  study  the  effect 
of   parameter  variations.      These  models   are   showm    in   Fig. 3.1. 

For  instance,  when  one  apply  the  theory  for  one 
parameter  of  interest,  say  A1,  the  sensitivity  equations  can 
be  obtained   from 


At 


rv-\i 


A. 


L    J 


l«i 


■\t 


[1 

r    ^ 

*z 

.  ^A 

y» 

+s; 

A3 

N 

xs 

L  .. 

M 

(3.  16) 


Similar    procedure  can   be   done      for   the  other  parame- 
ters   A2,    A3,   A4. 
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Figure  3.1    Nominal  and  Sensitivity  Models,  / 
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C.  SENSITIVITY    ANALYSIS 

A  computer  program  using  [Ref.  4]  was  written  in  order 
to  simulate  simultaneously  the  nominal  and  the  sensitivity 
models  as  shown  in  the  appendix  D.  A  step  was  applied  as 
input   of  the  nominal   system. 

The  number  of  equations  solved  are  6  for  the  nominal 
model  and  24  for  the  sensitivity  model.  Each  parameter  was 
varied  simultaneously   10%  from  all  the   nominal  values. 

D.  ANALYSIS  OF    -PARAMETER  VARIATIONS 

The  results     are   plotted  in  Fig. 3.2     to   3.6   and  Table  I 

Each    state     variable   output    is      plotted   separately  with  the 

correspondent  four  output  sensitivity  functions.  By  means  of 
these   plots, the    following  observations   can   be  made: 

The   plot  of  a,  :    =  *sTa "*"U-     >    j=  1,2,3,      and   *    in    fig.  3. 2 

indicates     that    a      parameter  change      A  Aj0( j=1,2, 3,      and     4) 

primarily  affect  the   rise  time   and  overshoot    of  X1(t)      since 

,111*     Z12,     ,113,     ,114,     are    largest      at   the  time    where   these 

effects   in   X1   occur. 

o  x  Z 
The   plot  of  /l^s   "^Sa"~'a*    '      j=1,2,3,      and   %    in    Fig. 3. 3 

indicates  that  X2(t)  is  little  affected  due  to  A  Alo  and 
fc>  A20  ana"  strongly  affected  in  the  rise  time  due  to  param- 
eter variations  ^A3o  and  ^&^.  The  overshoot  and  steady 
state  are  little  affected  due  to  parameter  variations 
A  A j*    (j=1,2,3,    and   4)  . 

The   plot  of  ,1,.=   ■55rr*-U*    *      j=1,2,3,      and  '4    in   Fig.3.4 
34  l)Aj       *\)0 

indicates  that  X3(t)  is  little  affected  in  +he  rise  time  due 
to  parameter  variations  A  ( j=1 ,2, 3, and  4).  The  overshoot 
and  steady  state  are  little  affected  due  to  parameter  varia- 
tions A  A,  0      and  AK2o  and   strongly     affected   due   to    h  A  30  and 

AA3T>- 
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The  plot  of  4.4;  ■  ^"/T-  'a'o'  j=1»2r3,  and  4  in  Fig.  3. 5 
indicates  that  X4(t)  is  little  affected  in  the  rise  time  and 
overshoot  due  to  parameter  variations  AAlo  and  ^  A  zo  and 
strongly  affected  due  to    ^  A3o  and    ^  A4o* 

The  plot  of  Z5j  =  s«~?—  K.  ,  3=1,2,3,  and  4  in  Fig.3.6 
indicates  that  X5(t)  is  little  affected  in  the  rise  time  due 
to  parameter  variations  A  \^Q  (js1,2,3,  and  4).  The  over- 
shoot is  little  affected  due  to  A  A  and  A  a2o  and  strongly 
affected  due  to    ^A3t)  and  A  A^0. 

The  plot  of  X(,j-  tT/^'W  '  3s1'2'3'  and  "  in  Fi<?-3-7 
indicates  that  X6(t)  is  little  affected  in  the  rise  time  due 
to  parameter  variations  £*  Aj  0  (j=1#2,3,  and  4)  .  The  over- 
shoot and  steady  state  are  little  affected  due  to  A  A(0  and 
4  A 20  and  strongly   affected   due  to  A  a3o  and   A  A^0. 

Table  I  shows  the  above  analysis  in  a  condensed  way  to 
give  a  general  picture  of  all  statss  and  sensitivity  func- 
tions  with   the   correspondent  effect   as  a   function   of   time. 

E.       PARAMETER-INDUCED  OUTPUT    ANALYSIS 

If     Zi«<   «    x0       ,the     actual      output   can      be    written     as 
(Eqn.2.29)  . 


y(t,*)    -  y<t,o^)     ♦    CTft,^)^^  (3.17) 


The  computer  program  as  shown  in  the  Appendix  D  gives 
the  actual  output  when  10%, 30%  and  H0%  of  variation  from  the 
nominal     value      of      each   parameter     is     assumed.  Figs. 3. 8 

through  3.13  show  the  actual  output  when  each  parameter  is 
varied   10%    from   the    nominal    value. 

Fig.  3. 8  indicates  that  the  overshoot  and  rise  time  of 
the  actual  and  nominal  output  of  X1  are  strongly  affected 
which   is  in  agreement  with    the   previous    analysis. 
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Fig. 3. 9  indicates  that  the  steady  state  and  overshoot  of 
the  actual  and  nominal  output  of  X2  are  primarily  affected 
which   is  in  agreement  with    the   previous   analysis. 

Fig. 3. 10  indicates  that  the  steady  state  and  overshoot 
of  the  actual  and  nominal  output  of  X3  are  primarily 
affected   which    is  in  agreement   with  the   previous    analysis. 

Fig. 3. 11  indicates  that  the  overshoot  of  tha  actual  and 
nominal  output  of  X4  is  primarily  affected  which  is  in 
agreement   with   the    previous   analysis. 

Fig. 3. 12  indicates  that  overshoot  of  the  actual  and 
nominal  output  of  X5  are  primarily  affected  which  is  in 
agreement   with   the    previous    analysis. 

Fig. 3. 13  indicates  that  the  steady  state  and  overshoot 
of  the  actual  and  nominal  output  of  X6  are  primarily 
affected  , having  little  effect  in  the  rise  time  which  is  in 
agreement   with   the   previous    analysis. 

Figs. 3. 14  through  3.19  and  Figs. 3. 20  through  3.25  show 
respectively  the  actual  output  for  30%  and  '4  0%  of  the 
nominal   value. 

From  the  plots  one  can  see  that  for  small  parameter 
variations  the  parameter-induced  output  error  is  negligible 
and  when  the  variation  becomes  large  as  30%  or  $0%  one  note 
that  the  error  becomes  pronounced  and  that  modeling  is 
starting  to  break  down.  This  agrees  with  the  assumption 
made    in   the  derivation  of  the   Eqn.2.28. 

F.       COUPLED    ROLL-TAI    AUTOPILOT    AN ALTS IS 

1-      Linear   Equations  of   the  Nominal   System 

From  the  block  diagram  of  the  Fig.B.1  in  the 
appendix  B  one  can  get  the  following  state  equations  of  the 
coupled  roll-yaw  autopilot. 


U6 


X1  ■  B  Conv  (A4  X3  +  A3  X8  ♦  A5  X11)  (3.18) 


X2  =  C  Conv  (A7  X3  ♦  A8  X8  +  A6  X11)  (3-19) 


X3  =  -  X1  -  (ALPHAB/Conv)  X2  (3.20) 

♦  KB  A  A2  X3  ♦  KB  A  A1  X11 
X4  =  -  8  X4  -  17.6  X12  ♦  17.6  PHC  (3.21) 


X5  =  -  (0.755/Conv)  X2  -  C  D  A7  X3  (3.22) 

♦     (0.755   -    8    D)     X4   -    5   X5   -    C    D   A8  X8    -    C    D   A5    X1  1 
17.6    X12   +17.6    PHC 

X6    =   C   E    A7   X3   -    6    X6    +   E    C    A8    X8    ♦  C    E    A6    X11  (3.23) 


X7    =   -   F    KC    (0.755/Conv)     X2  (3.24) 

-  (D    ♦    E)     F    KC   C    A7   X3    +    P    KC    (0.755   -    8    D)     X4 
+    KC(15   -    5   F)     X5    ♦    KC(6F    -15)    X6    -     15    X7 

-  (D    ♦    E)     F   KC    C   A8    X8    -    (D    ♦    E)     F    KC   C    A6    X11 

-  F    KC    D    17.6    X12    ♦    F    KC   D    17.6    PHC 

X8    =    188.4  Conv   X7    -    188.4    X8  (3.25) 
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19  ■  (K1  A  A2/t1)  X1  -  X9  /  T1  +  (K1  A  A1/T1)  X1 1  (3.26) 


X10  =-  (K2/Conv)  X1  -  (K2  H  ALPHAB/Conv  X2 


(3.27) 


U2/10)  (K1    A    A2/T1)    ♦    B    44    -   H    ALPHAB    C    A7    X3 
♦     (K2/10) (B  A3    -   H    ALPHAB    C   A8    X8    ♦    K2(1 

-  1/(10    T1))     X9    *     (K2/10)  (  (K1    A    A1>/T1    +    B    A5 

-  H    ALPHAB   C    A6)    X1 1 


X11    =    188. U   Conv    X10   -    188. 4    X1 1 


(3.28) 


X12   =    X2    /  Conv 


(3.29) 


Eqn.3.30      gives     the  matrix      representation     of      the 
state  variables    above  mentioned. 
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(3.30) 
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The   correspondence    of   the   states   is: 
X1   =   r,X2    =   p,X3   =   (i     , 
X4    =    X    ,    X5    =    11     ,    X6    =    X1    , 
X7    =   Sft.c     ,    X8   =    S*.  ,    X9    =    Y    , 
X10    =  &Yc     '     X11    =  SY    ,    X12    =   <b 

The  definition  of  the  constants  C1  throagh  C25  are 
given  in  appendix  B.  The  parameters  of  interest  for  this 
system   are   given   by    A1   through   A8    which   are    : 

A1    =   CYr      ,    A2    =    Cv        ,    A3    =   CH,       ,    A4    =    C         ,     A5    =    C  , 

A6    =   C0  ,    A7    =    C»  ,    18    =   C,       . 

<SY  aP  **«, 

2.      Sensitivity    Equations 

As  showed  in  the  previous  analysis  hare  one  can 
apply  the  same  procedure  using  tha  TRAJECTORY  SENSITIVITY 
EQUATION   as   given   in   Equ.2. 13. 


^i     =    Ao\i*    V Lxo  *    V L*(t)     »A;(^    =0     (3.31) 


For  the  purpose  of  this  procedure  again  one  can  see 
in   Egn.3.30   the   correspondent   matrices   A0  and  BQ.      _ 

The  partial  derivatives  ^-x-.*"|-  and  the  £rTr-\~~  are 
evaluated  considering  the  parameter  A  of  interest.  In  this 
case  they  are  A1 ,A2 , ... , A8,  that  are  respectively  the  aero- 
dynamic coefficients  present  in  the  coupled  roll-yaw  auto- 
pilot. 

Once  again, applying  the  partial  derivatives  with 
respect  to  the  parameters  one  can  obtain  eight  matrices 
respectively  as  found  , similarly,  in  the  previDus  case  of 
the    uncoupled  pitch    autopilot.         __ 

Here, one  can  see  that  x- — -  =  0  because  B  is  inde- 
pendent  of   the    parameters  of  interest. 

In   terms  of    components  one  obtains 
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(3.32) 


i=1,2, 12 

j=1,2, 8 

Here,one  can    see   that   eight      models  are  necessary   to 

study  the  effect   of    parameter     variations.      These    models   are 

shown    in  Fig. 3. 27.       For  one      parameter  of   interest,      say   A1, 

one    have  the  following  sensitivity  aquations    : 
!"•       1    r        -ir  ir         "1    r  .  1 

\i  •  * 


*3f 
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« 

*81 
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111 


At 
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*2 

x< 
*5 

X7 

XB 
Xg 

*11 


(3.33) 


Similar    procedure  can   be   maie    for   the   others   parame- 
ters   A2   through    A8. 


G.       SENSITIVITY    ANALYSIS 

In  order  to  simulate  simultaneouly  the  nominal  and  the 
sensitivity  models, a  computer  program  was  written  as  shown 
in  the  appendix    5. 

For  analysis  purpose,each  parameter  was  varied  simulta- 
neously   10%  from  the    nominal  value. 
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1-      Analysis   of  <X -Parameter   Variations 

The      plots    of       A,.'    =     -^T"1""~Ia-     *         j=1,2,...,8      in 

J  0  <\i       Ai  o 

Fig.  3.27  and  3.28  indicate  that  a  parameter  change  ^A-  , 
(j*1  ,2,.  •• ,  8)    primarily  affect   the   overshoot   of   JC1(t). 

The  plots  of  A-2j'  =  ^"'ajo'  j«1,2,...,8  in 
Fig.  3.29  and  3.30  indicate  that  a  parameter  change  £Ajor 
(j=1 ,  2,. .. ,  8)    primarily  affect   the   overshoot   of   X2(t). 

The  plots  of  ^Bj  =  ^!^--lAjo'  j=1,2,...,8  in 
Fig. 3.31  and  3.32  indicate  that  a  parameter  change  AAjor 
(j  =  1 ,  2r. .. ,  8)    primarily  affect  the  overshoot   of   X3(t). 

The  plots  of  A^:  =  m\j?~~  'aT  '  j=1r2r...,8  in 
Fig.  3.33  and  3.34  indicate  that  a  parameter  change  ^AjQ, 
(ja1  »2#. ..  ,8)    primarily   affect   the   overshoot    of   X4(t). 

The     plots   of     A5j     =     ^gXff.-i  j=1,2,...,8     in 

Fig.  3.35  and  3.36  indicate  that  a  parameter  change  ^Aj0, 
(j=1  #2,. .. ,  8)    primarily  affect   +  he   overshoot   of   X5(t). 

The  plots  of  \^  =  ^a^'Ajd'  j=1,2,...,8  in 
Fig.  3.37  and  3.38  indicate  that  a  parameter  change  AA^( 
(j-1 9  2,. ..  r  8)    primarily  affect   the  overshoot   of   X6(t). 

The  plots  of  A  7j*  ■  tsr  TV"- 1*-  #  *  j=1»2,...r8  in 
Fig.  3.39  and     3.40    indicate      that   a      parameter  change  AAj0, 


(j  =  1  ,  2r. .. ,  8)    primarily  affect^_the   overshoot   of   X7(t). 

Aj"       'A3  o 

Fig.  3.41   and     3.42    indicate      that  a      parameter  change  ^Ay0, 

(j*1#2f  ••  ,8)    primarily   affect  the  overshoot    of   X8(t). 

The     plots   of  A«i  •'      ■    ^-r^'-lA*    »        j-1,2,...,8     in 

J  vJ  o  Aj         Aj  o 

Fig.  3.43  and     3.44   indicate      that   a     parameter  change  ^AJQ, 

(j-1  #2r. .. ,  8)      primarily  affect  the      overshoot   and   rise  time 

of  X9(t)  .  — N 

v        .  (JXio 

The      plots    ofAioj      =    "J™- — 'AT  '         j=1,2r...,8      in 

Fig. 3.45  and     3.46    indicate      that  a      parameter  change  A  A,    , 
(j=1  r  2r. .. ,  8)    primarily  affect  the   rise   time   of   X10(t). 


The     plots   of    Aa^'     =    -&fz 
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v  U  X  i| 

The     plots   of    AUj     ■   — -----|    .     *         js1»2, ,8      in 

Fig.  3.47  and  3.48  indicate  that  a  parameter  change^  A3  0, 
(j=1 ,2,. ..  ,8)    primarily   affect  the  overshoot   of   X11(t). 

The  plots  of  Ajz^j*  =  ^~£-  U^o*  j=1r2,...r8  in 
Fig. 3.49  and  3.50  indicate  that  a  parameter  change  A Aj0  , 
(j*1 ,2,.  ••  ,8)    primarily  affect   the  overshoot    of   X12(t). 

Table  II  and  III  show  the  above  analysis  in  a 
condensed  way  to  give  a  general  picture  of  all  states  and 
output  sensitivity  functions  with  the  correspondent  effect 
as   function  of   time. 

2-      Parameter- Induced  Output    Ani.i.Z§i.s_ 

As  Shown  previously,  if  ^  oc  <<  oc0  ,  the  actual 
output  can    be   written  as 


f  <t,K)    t  y(t,og     ♦<T(t,og    A<x  (3.34) 


The  computer  program  given  in  appendix  E  was  written 
for  simulating  the  system  when  10%,  30%,  and  43%  of  varia- 
tion from  the  nominal  value  of  each  parameter  is 
assumed. Fig. 3. 51  through  3.62  give  the  actual  output  when 
each    parameter    is  varied   10%  from  the   nominal   value. 

Fig. 3. 63  and  3.64  give  the  actual  output  for  30%  of 
variation  from  the  ncminal  value  of  aach  parameter  assumed. 
Fig. 3.65  and  3.66  give  the  actual  output  for  40%  of  varia- 
tion from  the  nominal  value  of  each  parameter  assumed. 
These  plots  show  the  output  of  the  state  variables  X3  and 
X11  that  present  strong  deviations  just  to  give  the  behavior 
of  the  system  when  parameter  variations  are  not  small.  One 
notes  that   modeling    is  starting  to   break   down. 
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Figure   3.2        Sensitivity  of    X1    with   Respect   to  A1,A2,A3,A4. 
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Figure  3.3        Sensitivity  of    X2   with   Respect   to  A1,A2,A3,A4. 
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Pigure  3,4        Sensitivity  of    X3   with   Respect  to  A1,A2,A3,A4. 
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Figare  3.5        Sensitivity  of    X4   with   Respect  to  A1,A2,A3,A4. 
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Figure  3.6        Sensitivity  of    X5   with   Respect  to  A1,A2,A3,A4. 
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Figure   3.7       Sensitivity  of    X6   with   Respect  to  A1,A2,A3,A4. 
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Figure  3,8        actual    and   Nominal  Output  of   X1     (10X    variation) 
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Figure  3.9        ictual   and  Noainal  Output  of   X2    (10X   variation) 
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Pigure  3.10   Actual  and  Nominal  Output  of  X3  (10%  variation) 
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Figure  3.11   Actual  and  Nominal  Output  of  X4  (10 X  variation) 
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Figure  3.12   Actual  and  Nominal  Output  of  X5  (1Q%  variation) 
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Pigure  3.13       Actual   and  Nominal  Output   of   X6    (13%   variation) 
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Figure  3.14   Actual  and  Nominal  Output  of  X1  (30%  variation) 
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Pigure  3.15   Actual  and  Nominal  Output  of  X2  (33%  variation) • 
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Figure  3.16        Actual   and  Nominal  Oat  put   of   X3    (30%  variation) 
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Pigure  3.17       Actual   and  Nominal  Output   of  14    (33%   variation) 
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Figure  3.18        Actual   and  Noninal  Output   of   X5    {30%  variation) 
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Figure  3.19       Actual   and  Nominal  Output   of   X6    (33%   variation) 
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Figure  3.20   Actual  and  Noainal  Output  of  X1  (43%  variation) 
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Figure  3.21   Actual  and  Nominal  Output  of  X2  (15%  variation) 
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Pigure   3.22       Actual   and  Nominal  Output   of   X3    (43  36   variation) 
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Figure   3.23        Actual   and  Nominal  Output   of   X4    (43%  variation)  . 
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Figure  3.24        Actual   and  Nominal  Oat pat   of   X5    (40 X  variation) • 
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Figure  3.25        Actual   and  Nominal  Output  of   X6    (13%  variation) 
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Figure  3.26    Roll-Yaw  Nominal  and  Sensitivity  Hodels- 
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Figure   3.27        Sensitivity  of  X1   with    Respect   to   &1rA2rA3,A4, 
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Figure  3.28   Sensitivity  of  X1  with  Respect  to  *5,A6,A7,A8. 
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Figure   3.29        Sensitivity  of  X2  with    Respect   to   MrA2,A3#A4. 
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Pigure  3.30   Sensitivity  of  X2  with  Respect  to  i5,&6,A7rA8. 
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Figure   3.31        Sensitivity  of  X3   with    Respect   to   M,A2rA3,A4. 
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Figure  3.32   Sensitivity  of  X3  with  Hespect  to  &5,A6rA7,A8. 
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Pigure  3.33        Sensitivity  of  14  with.    Respect   to   M,A2,A3rA4. 
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Figure  3.34   Sensitivity  of  X4  with  Respect  to  &5,A6,A7,A8. 
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Figure  3,35   Sensitivity  of  X5  with  Respect  to  M,A2,A3,A4. 
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Pigure  3.36   Sensitivity  of  15  with  Respect  to  &5#A6,A7, A8. 
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Figure  3,37   Sensitivity  of  X6  with  Respect  to  M,A2,A3,A4. 


89 


o 

z 


Q      G      <       +      X 


i/)       to       r»       co 
to      <o       co       to       to 

X        -1        -J        -»        -1 


Saa 


02  •'o  gi  -'o  or'a  scro  ooJo        so-d-        ai-d-        srd-        02-d- 


90 -a  ho'o 


S91 


oco  02  -q  oca 


00 '0  01-0-  02'0- 

991 


OE'O-  Oh'O-  OS'O- 


EO   0  10*0  00-0-  20'0- 

LQ1 


E0*0-  WO-  90*0- 


oo  -I  08-'o  09  "'o  ovo  02 -'a 


00-0  02 -0- 

991 


Oft'O-  09*0- 


Figure  3.38   Sensitivity  of  X6  with  Respect  to  &5,A6,A7,A8. 
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Figure  3,39   Sensitivity  of  X7  with  Respect  to  \ 1,A2, A3, A«*. 
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Figure  3.41        Sensitivity  of  X8  with   Respect   to  &1r&2,A3rA4. 
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Figure  3,42   Sensitivity  of  X8  with  Respect  to  &5,A6,A7,A8. 
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Figure  3.43       Sensitivity  of  X9   with    Respect   to   M,A2,A3rA4. 
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Figure  3.44   Sensitivity  of  X9  with  Respect  to  K5, A6r A7, A8, 
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Figure  3.U5        Sensitivity  of  X10    wita   Respect   to   A1,A2,A3,AU 
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Figure  3.46   Sensitivity  of  X10  with  Respect  to  A5,A6,A7,A8 
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Figure  3.47        Sensitivity  of  X11    with   Respect   to   A1,A2VA3,A4 
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Figure  3.48   Sensitivity  of  X11  with  Respect  to  A5,A6,A7,A8 
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Pigure  3.49        Sensitivity  of   X12   with   Respect   to    A1,A2,A3,AU 
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Pigure  3.50        Sensitivity  of  X12   with   Respect   to    A5,A6,A7,A8 
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Figure  3.51   Actual  and  Nominal  Output  of  X1  (10*  variation) 
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Figure  3.52        Actual   and  Nominal  Output   of   X2(10£    variation). 
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Figure  3.53        Actual  and  Nominal  Output   of   X3(10X    variation) 
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Figure  3-5'*        Actual   and  Nominal  Output   of   X4(10£    variation) 
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Figure  3,55        Actual   and  Nominal  Output   of   X5(10S    variation). 
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Figure  3.56        Actual   and  Nominal  Output   of   X6(10R    variation) 
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Figure  3.57        Actual   and  Nominal  Output   of   X7(10X    variation) 
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Figure  3.58   actual  and  Nominal  Output  of  X8 (10X  variation) 
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Figure  3.59   Actual  and  Nominal  Output  of  X9  (10X  variation) 
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Figure  3.60       Actual  and  Nominal  Output   of   X10(13%  variat 


ion) 
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Figure  3,61   Actual  and  Nominal  Output  of  X1 1 (13 %  variation). 
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Figure  3.62       Actual  and  Nominal  Output   of   X12(10%  variation) 
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Figure  3.63   Actual  and  Nominal  Output  of  X3  (30X  variation) 
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Figure  3.64        Actual   and  Noainal  Output   of  X11(39%  variation) 
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Figure  3.65   Actual  and  Nominal  Output  of  X3  (40*  variation) 


119 


>-  >- 


OS'I  00*1  OS'O  OO'O 


Saa 


OS't  OO't 


OS'O-  OO'I- 

II JL 


os''o         otr'o 


OS'O-  OO'I- 

StU 


~n 


OS'l"  00'2-  OS'2- 


OS-I-  00'8-  0S*2- 


Figure  3,66        Actual   and  Nominal  Output   of   X11(U3*   variation) 
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IV.  APPLICATION  OF  SENSITIVITY  ANALYSIS  TO  NONLINEAR  SYSTEMS 

A.  INTRODUCTION 

The  three  dimensional  nonlinear  system  used  here  for 
sensitivity  analysis  purposes  is  presented  in  appendix  C. 
The  missile  configuration,  size  and  mass  properties  are  the 
same  as  those  used  in  the  linear  system  and  were  presented 
in   appendix   A. 

In  section  B  the  sensitivity  equations  with  respect  to 
parameters  of  nominal  system  is  shown  and  in  section  C  the 
sensitivity  equations  with  respect  to  parameters  that  showed 
to   be   of  more  effect    in  the    previous   analysis   is    presented. 

Section  D  shows  the  results  of  the  trajectory  sensi- 
tivity equations  when  step  inputs  are  applied  at  specific 
time    as   mentioned  in   appendix   C. 

Section  E  gives  the  parameter-induced  outpit  analysis 
when    each   parameter    is  varied    10%   from   the   nominal   value. 

B.  NONLINEAR  EQUATIONS  OF  THE  NOMINAL  SYSTEM 

From  the  block  diagram  of  Fig. C.1  and  Eqn.  C.3  through 
C.21  as  given  in  appendix  C  one  have  the  following  nominal 
equations. 

• 

It    *   -  C4  It    -    C02    C4    A2  (4.1) 


X2  -  -  C5  X2  -  C6  NZC  ♦  C6  C19  C03(X7)  ♦  C7  X1       (4.2) 


121 


X3  =     (C5    C8    -    C9)     X2    ♦    C6    C8    NZC  (4.3) 

-  C6    C8    C19    Cos(X7)    -   C7    C8    X1    ♦  (C8/Conv2)     X1 7    X18 

♦  C1    C8    A1    ♦     (C9/Conv)     X5 

X4  =   -    C3    X4    ♦    C3    Conv    X3  (4.U) 


X5    =    (X17   X18)/Conv   +    C1     Conv    A1  (4.5) 


X6    =   X5   -    KB    C02    A2  -    (X16    X18)/Conv  (4.6) 


X7    =   X5   Cos  (X1  9) /Conv   -   X17    Sin  (X1  9) /Conv  (U.7) 


X8  -    -    C10   X8    ♦    (C12   -    C1  1    C13)     X10  (4.8) 

♦  C11    C14   PHC    -    C11    C14    X19    -    C11    CO  (A8    X16  ♦    A7    X15 

+  A3   X13)    -     (C12/Conv)    X18 

X9  =   -    X9/T1    ♦    K1/T1    C02(AA    X16    ♦    A6    X15)  (4.9) 


X10    =    -    C13    X10    +    C14    PHC    -    C14    X1 9  (4.10) 


X11    -    -   C15    X11    ♦    C16    C0(A8    X16    ♦    A7    X15    +    A3    K13)   (4.11) 
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X12    =    -    C17    X12    ♦    K(C17    -     (C10/C13)     X8  (4.12) 

(K/C18)  (C12    -    C11    C13)     X10    ♦    K    C1  1  (C14/C18)     P  HC 

-    K    C11     (C14/C18)     X19    -     (K/C18)     CO  A8(C11    -    (K/C18) 

CO    A3(C11    *■    C16)     X13    -    (K/C12)  (C1 8  Conv)     X18 
♦    K(C15/C18)     -   C17)    X11 

X13    =    -    C3   X13    ♦    C3   Conv    X12  (4.13) 


X14    ■     (K2/10)  (-    (X9/T1)     C02(AA    X16    ♦    A6    X15)  (4.  14) 

-  X15    (X18/Conv2)     ♦   C01(A9    X16    +    A4    X15    ♦    A5    K13) 

-  (KYP/Conv2)  (X5    -  KB   C02    A2   -    X16  (X1 8/Con v)    X18 

-  (KYP/Conv)     XX6    C0(a8   X16    ♦    A7    X15    ♦    A3    X13) 
♦     (K2/Conv)     X17    -    K2(KYP/Conv2)     XX6    X18 

X15   =    -   C3   X15   ♦    C3   Conv    X14  (4,15) 


X16   -    KB    C02(AA   X16   ♦    A6    X15)     +    X5(X18/Conv)    -    X17(4.16) 


X17    =   -    (X5   X18)/Conv   ♦   C01    Conv(&9    X16  (4.17) 

+    A4    X15    +    A5    X13) 
X18   =    CO    Coav(A8    X16   +    A7    X15    A3    X13)  (4.18) 


X19    ■    X18/Conv  (4. 19) 
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The  correspondence  of  nonlinear  state  vectors1  as  given 
in  the  above  equations  with  relation  to  the  nonlinear  system 
presented   in  Fig.C.1    in   appendix   C   are    : 

X1    =    X  r  X2  =    Y      ,  X3    =  <$f»c  ,  X4    =     Sp  ,      X5    =   q    , 

X6    =  <X  ,  XI  -  O      ,  X8   =    Y1    ,  X9    =    Y2    ,    X1  0    =    X1    , 

X11    =    X2  ,  X12  =  &RC,  X13   =    SR  ,  X14    =   Syc,    X1  5    =    8r  , 

X16   =  (*>  ,  X17  =   r      ,  X18   =   p      ,  and    X19   =  <)>  . 

The  parameters  of  interest  for  this  nonlinear  system  are 
given   by 

A1    =CU,$>),    A2    =C<V/»),    A3    =(LU),    A  4   =  Cfll*J  ,    aad   45    =^*). 

Definition  of  the  constants  CO  through  C19  are  given  in 
appendix  C. 

As  previously,  this  nomenclature  is  used  here  to  easily 
apply   the  sensitivity  theory. 

C.       NONLINEAR    SENSITIVITY   EQUATIONS 

From  the  sensitivity  theory  given  in  chapter  2  one  knows 
from    Egn.2.12   that    for  nonlinear   systems  one   has 


"Dx         ^H      ^x  n>{  ^X0 

.J— +    —J—    t    =    o  (4.20) 


^^      "c)x     T)*j'       D^'       1> 


<*J 


Here 


X    -    f  {X1,X2,  ..  .,X19,t,NZC,A1,A2,..  .  ,A5}  (4.21) 


lThe  state  vectors  and  parameters  in  the  noi linear  case 
have  no  correspondence  to  the  state  vectors  and  parameters 
used  in  the  linear  case. The  correspondence  can  be  found  in 
the  appendices  B  and  C. 
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From   Eqn.2.14   one   knows      that   the   trajectory   sensitivity 
vector   is  given    as 


A        =   --A—  |_A     ♦    ~ |     ,    A  (0)     =0  (4.22) 

OX       *0  J  Oo^     *0 


Applying  the  above  theory  one  obtains  tha  nonlinear 
sensitivity  equations  as  given  in  Eqn.4.24  through  S4e42, 
where     ,A        was    considered    the    parameter   to   be  varied. 


211    =    -   C4/111    -    C02   C4(DA26,161    ♦    DA24,141)  (4.23) 


,121    -    -   C5A21    -    C6  C19   Sin(X7)    ,171    ♦    C7,111  (4.24) 


,131    =    (C5   C8    -   C9)   ,121    ♦    C6    C8    C19    Sin  (X7)   ,171  (4.25) 

-    C7   C8  ,111    +   C8(x18/i171   +    (X17  A181)/(Conv2> 

♦    C1    C8(DA16/i61    ♦   DA  14  ,141    ♦    1.)     ♦    (C9  ,15 1)  /Conv 

,14  1    =   -  C3/141    ♦    C3  Conv  ,131  (4.26) 


,151    =     (X17,1181    ♦    X18/i17  1)/Conv  (4.27) 

+    C1    Conv(DA16  ,161    +    DA14,141    ♦    1.) 
/161    =,151    -    KB   C02(DA26/161    ♦    DA24,141)  (4.28) 

-     ((X 16/1181    ♦   X  18  .1161)  /Conv 
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,171    ■     (Cos(X19)  /151    -    X5    Sin  (X19)  ^1191)  /Conv 


(4.29) 


-     (Sin(X19)    ,1171    ♦   X17   Cos  (X19)  A.  1  71 )  /Conv 


,181    =    -    C10,181    ♦     (C12    -    C11    C13|   A101 


(4.30) 


-    C11    C14,1191    -    C11    CO     (A8/1161    ♦    A7,1151    +    \3^i131 
♦  /161(X16  Da8   ♦    X15   Da7    ♦   X13    DA13))     -    (C1 2  ,11  81)  /Conv 


,191    =    -,191/T1    ♦    K1/(T1    C02)   (AA^.161 


(4.31) 


♦    A6  /1151   ♦  ,161  (X16   DAA    ♦    X15    Da6)  ) 


,1101    =   -    C13,1101     -   C14/1191 


/ 


<<*.32) 


,1111    =   -    C15/1111    *   C16    C0(A8,1161    +    A7/1151 


(4.33) 


♦    A3/1131    ♦  /161(X16   DA8    +    X15    DA7    ♦    X13    DA3)  ) 


,1121    =   -    C17xL121     +   K(C17    -    C10/C18)   ,181 


<»-3U) 


+  (K/C18)  (C12    -    C11    C13)    /1101    -    K    C11     (C14/C18)   ^.191 

-  (K/C18)    C0(C11     ♦   C16)  (A8^.161    ♦    A7  ,1151    ♦    A3  ,1131 
+    A3  ,1131   «•  ,161  (X16   DA8    ♦    X15    DA7    ♦    X13    Da3)  ) 

-  K    C12/(C18    Conv)  ,1181    ♦    K    (C15/C18    -    C17),1111 


,1131    =   -    C3^131    ♦    C3    Conv^.121 


(4.35) 
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,1141    =    (K2/10)   (-    X9/T1    ♦     (K1/T1)     C02(AA  (4.36) 

,1161    +    A6/1151    ♦   a61(X16    DAA    ♦    X15    DA6)  )     -     (X15/1181 

♦  X18  Z151)  /(Conv2)    +    C01(A9Z161     +    A4,1151    ♦    A5^.131 

♦  ,161  (X16   DA9    ♦    X15    DA4     ♦    X13    DA5|  )     -    KYP/(Conv2) 
(X5  2.181    ♦    X18/151   -    KB    C02(A2/1181    ♦    X18(DA25^i61 

♦  DA24/141)     -(X182,1161    ♦    2    X18    X16/1181)     /Cotiv)) 
-     (KYP/Conv)  C0(    xx6/i161    ♦    X16,116  1)     ♦    A7,1161) 

+    A7(XX6/i151    ♦    X15/1161)    +    A3(XX6/i131    +    X13A161) 

♦  X16(XX6    DA8/161    +    A8,1161)     ♦    X15(XX6    DA7^161 

♦  A7  2.161)     ♦    X13(XX6    DA  3  ,161    ♦    A3  ,1161))) 

♦  (K2/Conv)    ,1171    +   K2/19  1   -    K2   KXP/(Conv2) 
(XX6  ,1181   ♦    X18  ,1161)) 


• 


,1151    =   -    C3  2.151    ♦    C3    Coiiv  2.141  (4.37) 


2.161    =   KB   C02(AA  ,161    ♦    A6   2151    +  2.61(  (4.38) 

X16    DAA    ♦   X15    DA6))    ♦    (X6/1181    ♦    X 18  ,16  1)  /Con?    -A171 
•1171    =   -     (X5  2.181    +  X18/151) /Conv  (4.39) 

♦  C01    Conv(A9,l161   ♦    A4/1151    +    A5  .2.131  +  ,161(X16    DA9 

♦  X15    DA4   ♦    X13    DA5)) 

,1181    =  CO   Conv  (A8  ,1161    ♦    A7,1151  (4.40) 

♦  A3  2.131    +  ,161(X16   DA8   ♦    X15    DA7    ♦    X13  DA3)) 

,1191   =/i181/Conv  (4.41) 
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The  same  procedure  above  must  be  made  foe  the  other 
parameters  that  most  affected  the  previous  analysis  and  will 
be  given  in  the  nonlinear  sensitivity  analysis  presented  in 
the    next  section. 

D.       NONLINEAR   SENSITIVITY   ANALYSIS 

A  computer  program  in  appendix  3  shows  the  simulation  of 
the  sensitivity  equations  with  respect  to  the  selected 
parameters. 

The  number  of  equations  solved  are  19  for  the  nominal 
equations  and  95  for  the  sensitivity  equations.  Here,  the 
parameters  were  selected  from  the  previous  two  cases  given 
in  Chapter  III  that  showed  to  have  most  effect  in  the  time 
response   of  the    systems. 

From     the     analysis     performed,  two      paraaeters      were 

selected  from  the  uncoupled  pitch  autopilot  and  three  from 
the  coupled  roll-yaw  autopilot,  which  has  been  shown  to  be 
more  sensible  to  parameter  variations.  These  parameters  are 
respectively  CN^  (A2)  and  Cm  (A4)  in  the  uncoupled  pitch 
autopilot  and  Q?  (A3),  Cn,  (A5)  ,  and  Cn  (\Q)  in  the 
coupled   roll-yaw   autopilot. 

In  order  to  compare  the  influence  of  the  parameter  of 
interest  in  the  nonlinear  case  some  state  variables  were 
chosen  to  be  analysed.      The    state   variables   selected   were: 

Prom   the   uncoupled  pitch   autopilot    : 

X3    *  Sp^  ,    X4    =    £p  ,    X5    =    q    ,    X6    =    X 

From   the  coupled    roll-yaw  autopilot    : 

X12   =  <SRc    ,     X13    =    &K     ,     X14    =    5yc    ,    X15    =     SY     , 

X16  =  (%       ,    X17    =   r   ,    X18  =   p    ,    and   X19    =    £     . 
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The  results  are  plotted  in  Fig- 4.1  through t  4.12  and 
Table  IV  and  V.  Each  state  variable  output  is  plotted 
separately  with  the  correspondent  five  output  sensitivity 
functions. 

From  the  plots,    the   following   observations  can   be   made: 

Fig.  4.1  indicates  that  the  state  variable  X3(6j>  )  is 
little  affected  in  the  rise  time  by  the  parameters  A1  and 
A2.  The  overshoot  is  strongly  affected  by  A5,  A1  and  A4  and 
little  affected  by  A2  and  A3.  The  steady  state  is  strongly 
affected  by   A1    and    little  affected  by   A2. 

Fig. 4. 2  indicates  that  the  state  variable  X4(5p)  is 
little  affected  in  the  rise  time  by  A1  and  A2.  The  over- 
shoot is  strongly  affected  by  A1,  A3,  and  A4  with  little 
effect  by  A2  and  A5 .  The  steady  state  is  strongly  affected 
by   A1    and  little   affected  by   A2. 

Fig. 4. 3  indicates  that  the  state  variable  X5  (  ^  )  is 
strongly  affected  in  the  rise  time  by  A1  and  A2.  The  over- 
shoot is  strongly  affected  by  A3r  A4,  A5,  and  little 
affected  by  A1  and  A2.  The  steady  state  is  not  affected  by 
parameter  variations. 

Fig. 4. 4  indicates  that  the  state  variable  X6  (  *<  )  is 
strongly  affected  in  the  rise  time  by  A1  and  A2.  The  over- 
shoot is  strongly  affected  by  A3,  k4,  and  A5.  The  steady 
state   is  little   affected  by    A1   and  A2. 

Fig. 4. 5  indicates  that  the  rise  time  and  steady  state  of 
the  state  variable  X12(6rcJ  are  not  affected  by  parameter 
variations.  The  overshoot  is  strongly  affected  by  A3  and  A5 
and    little  affected    by  A1,    A2,    and  A4. 

Fig. 4.6  indicates  that  the  rise  time  and  steady  state  of 
the  state  variable  X13(&j^)  are  not  affected  by  parameter 
variations.  The  overshoot  is  strongly  affected  by  A3,  A4, 
and    A5  and   little  affected    by   A1    and    A2. 

Fig. 4. 7  indicates  that  the  rise  time  and  the  steady 
state  of  the  state   variable    X14(Sy  )    are  not   affected   by  the 
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parameter  variations.  The  overshoot  is  strongly  affected  by 
A4   and   A5   and   little   affected   by    A1,    A2,    and    A3. 

Fig. 4. 8  indicates  that  the  rise  time  and  the  steady 
state  of  the  state  variable  X15(6y  >  are  not  affected  by  the 
parameter  variations.  The  overshoot  is  strongly  affected  by 
A3,    A4,    and   A5    and    little  affected   by   A1,    and   A2. 

Fig. 4. 9  indicates  that  the  rise  time  and  the  steady 
state  of  the  state  variable  X16(a  )  are  not  affected  by  the 
parameter  variations.  The  overshoot  is  strongly  affected  by 
A3r    A4,    and   A5    and    little  affected  by   A1,    and  A2. 

Fig. 4. 10  indicates  that  the  rise  time  and  the  steady 
state  of  the  state  variable  X17(  \-  )  are  not  affected  by  the 
parameter  variations.  The  overshoot  is  strongly  affected  by 
A3,    A4,    and  A5    and    little   affected   by   A1,    and   A2. 

Fig. 4. 11  indicates  that  the  rise  time  and  the  steady 
state  of  the  state  variable  X18(  f»  >  are  not  affected  by  the 
parameter  variations.  The  overshoot  is  strongly  affected  by 
A3,    A4,    and   A5    and    little   affected   by   A1,    and   A2. 

Fig. 4. 12  indicates  that  the  rise  time  and  the  steady 
state  of  the  state  variable  X19(  <f  |  are  not  affected  by  the 
parameter  variations.  The  overshoot  is  strongly  affected  by 
A3,    A4,    and   A5    and    little  affected   by   A1,    and   A2. 

E.       PABAMETEB-INDOCED  OUTPUT   ANALYSIS 

As  shown  previousy,  if  A<*«  o(Q  r  the  actual  output  can 
be   written   as 


a 


♦  <r 


V     (t,*)     =    y(t,  o^)     ♦    U     (t,c*0  )    4<X  (4.42) 


The  computer  program  given  in  appendix  G  was  written  for 
simulating  the  system  when  variations  from  the  nominal  value 
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of  each  parameter  is  assumed.  Figs. 4. 13  throught  4.24  giva 
the  plots  of  the  actual  and  nominal  output  of  the  state 
variables  that  are  present  in  the  unclouped  pitch  autopilot 
and  coupled  roll-yaw  autopilot  when  10%  of  parameter  varia- 
tion is  assumed.  The  others  state  variables  are  not  given 
here    because  they  are   just    intermediate   state   variables. 

these  results  confirm  the  theory  presented  in  chapter  2 
when  small  parameter  variations  are  assumed.  Figs. 4. 25  and 
4.26  present  the  actual  and  nominal  output  of  the  state 
variables  X3  and  XI 1  that  showed  pronounced  variations  when 
30R  of  parameter  variation  is  assumed.  Figs.  4.27  and  4.28 
present  the  actual  and  nominal  output  of  X3  and  X 1 1  when  40% 
of  parameter     variation  is   assumed.  From   these      polts   one 

note   that  the  modeling  is  starting  to   break   down. 
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Pigure  4.1        Sensitivity  of    13  with   respect   to  A1 , A2r A3,A4, A5 
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Pigure  4.2       Sensitivity  of    XU   with   respect   to  A1 , A2,A3,A4r A5 
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Figure  4,3       Sensitivity  of    X5   with   respect   to  A1 ,A2,A3,A4,A5 
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Figure  4.U 


Sensitivity  of    16   with   respect   to   A1 ,A2, A3,A4, A5 
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Figure  4.5 


Sensitivity  of    11 2  with    respect   to   k 1 , \2,  A3, A4, A5, 
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Figure  4.6        Sensitivity  of    X13   with    respect   to   A  1, A2, A3r A4r A5. 

137 


□      0 

4     + 

X 

O 

a 
z 

III 

J    J    =» 

zf          ID 

• 

a 
r*rj 

O 

8    w 

- 

1 

3 

1 

l 

9 

J 

^ 

\ 

7* 

~» — " 

—zz^z 

Ss 

1— ■■ 

> 

<L 

~~^»^s=,,?e 

3 

'a 

cz 

<^. 

tl-n  ■  „,_ 

<r 

—» •" 

—              >*- 

3 
*«» 

1U 

« 

i 

i 

Id 

X 

i 

i 

> 

i 
i 

i 

i 
i 
i 
i 

3 
'•» 

8 

"ni 

8 

3 

90J0              *J0 

wo 

20  r0 

00 ! 

0            20'( 

TilX 

-           10'C 

00  "0- 

I0*( 

P«?H 

. 

SOJ0              iVo 

so  a 

w'o 

so- 

0              CO' 

inn 

0          OO'C 

V       zo-d- 

so-d- 

sr'o          zi'o 

01    0 

io-'o 

so- 

0             20' 

2fin 

0        oo-c 

-        eo-d- 

so-o- 

os-'z         oo-'z 

Os'l 

00 -1! 

os- 

0            00' 

etin 

o        ace 

BO'f- 

QS-T'- 

os't         w'l 

os-'o 

oo-a 

os  -c 

-        oo't 
nun 

-       os-r 

-       oo-z- 

DB'i 

_ 

00 -'8              OBJi 

CO '9 

osn 

GO- 

S          os- 

STin 

I         oo- 

0          U'1- 

00*2- 

• 



Figure  4.7       Sensitivity  of    X14  with    respect   to   M,  A2,  A3,  A4,  A5. 
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Pigure  4.8        Sensitivity  of    X15   with    respect   to   kl, A2, A3, A4, A5. 
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Figure  4.9        Sensitivity  of    X16   with    respect   to   M, A2, A3, A4, A5. 
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Figure  4.10       Sensitivity  of  X17   with   respect   to    A1, A2, A3rA4, A5, 
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Figure  a. 11        Sensitivity  of  X18    with   respect   to    A1,  A2, A3,A4r  A5 
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Figure  4.12        Sensitivity  of  X19    with   respect   to    A1, A2, A3rA4, A5. 
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Figure  4,13   Actual  and  Nominal  Output  of  X3  (13%  variation) 
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Figure  4.14   Actual  and  Nominal  Output  of  X4  (19%  variation) . 
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Pigure  1.15       Actual   ana  nominal  Output   of   X5    (13%  variation) 
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Figure  4.16   Actual  and  Nominal  Output  of  X6  (13%  variation) 

147 


B      Q 

LEGEND 
Y12 
Y12S 

a 
a 

1 1 

_o 

9*.  00 

Sec 

a 
a 

a 

o 

5="" 

a 

a 

o 

"(O 

a 

a-= 

"*~+ 

i 

"i/i 

UJ 

3C 

'^ 

a 
a 

"cri 

o 
a 

"ni 

a 
a 

a 
9 

oe  •  a          sz  •  o 

02'0 

9ro 

IfO 

90 '0 

2IJL 

20  '0            £D-0- 

BO'C 

p»h 

oe  •fa         sz  -'o 

020 

9['0 

Il'O 

90*0 

S2U 

20-Q           EQ-O- 

80*0 

.    _                                — 

■             

Pigure  4.17   Actual  and  Noninal  Output  of  XI 2  (10%  variation) 
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Figure  4.18        Actual   and  Nominal  Output   of   X13    (10%   variation) 
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Figure  4.19   Actual  and  Noainal  Output  of  X14  (10%  variation). 
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Figure  4,20        Actual   and  Nominal  Output   of   X15    (10%  variation) 
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Figare  4-21   Actual  and  nominal  Output  of  X16  (10%  variation) 
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Figure  4,22       Actual  and  Nominal  Output   of   117    (10%  variation) 
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Lqure  4.23       Actual  and  No.inal  Output   of  X18    (10%  variation) 
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Figure  4.24        Actual   and  Nominal  Output   of   X19    (10%  variation) 
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Figure  4.25   Actual  and  Nominal  output  of  X6  (39%  variation) 
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Figure   U.26        Actual   and  Nominal  Output   of  X16    (30%   variation) 
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Figure  4.27        Actual   and  Nominal  Output   of   X6    <«0%   variation) 
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Figure  4.28        Actual  and  Nominal  Output  of   116    (10%  variation) 
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?,     CONCLUSIONS    AND   RECOMMENDATIONS 

A.       CONCLUSIONS 

The  results  of  the  nonlinear  3-D  parameter  sensitivity 
analysis  presented  in  Chapter  IV  verify  the  linear  analysis 
of  Chapter  III-  These  conclusions  can  be  quickly  made, 
using  Tables  I,  II,  III  and  17  that  give  a  brief  review  of 
the  time  response  of  the  linear  and  nonlinear  system. 
Another  means  of  comparison  of  the  linear  and  nonlinear 
analysis  is  to  use  the  figures  that  give  a  precise  visuali- 
zation  of   the  time   response. 

The  comparison  of  the  linear  and  nonlinear  analysis  can 
be   done   ty    means   of    the    plots  as   follows: 


B.       EFFECT    ON   STATE    VABIABLES    DOE   TO    PARAMETER   VARIATIONS 

1-      Effect      on      State      Variable        Due      to      Variation     in 

«n 

Figs. 3. 7  and  U.1  (  Ac^  and  A^j  )  show  that  ^f>c  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot   and  steady   state. 

Figs. 3. 4  and  4.2  (  X3.4  and  X4  j  )  show  that  6  f>  is 
little  affected  in  the  rise  time  and  strongly  affected  in 
the   overshoot   and  steady  state. 

Figs. 3. 2  and  4.3  (A^  and  A51)  show  that  g  is 
strongly  affected  in  the  rise  time,  little  affected  in  the 
overshoot   and  not  affected    in   the   steady  state. 

Figs. 3. 3  and  4.4  (  A24  an:i  A^l  )  show  that  «<  is 
strongly  affected  in  the  rise  time,  little  affected  in  the 
overshoot   and  steady   state. 


162 


2.  Effect      on      State      Variable        Due      to      Variation     in 

M 

Figs. 3. 7  and  4.1  (  Afo2  and  ^32)  show  that  &PC  is 
little  affected  in  the  rise  time,  overshoot  and  steady 
state. 

Figs.  3.4  and  4.2  (A32  an:i  ^42  )  show  that  Sp  is 
little  affected  in  the  rise  time,  overshoot  and  steady 
state. 

Figs. 3. 2  and  4.3  (  A12  and  /V52.)  show  that  q  is 
strongly  affected  in  the  rise  time  and  little  affected  in 
the   overshoot   and  steady  state. 

Figs. 3. 3  and  4.4  (^22  an:i  ^62)  show  that  ©<  is 
strongly  affected  in  the  rise  time  and  little  affected  in 
the   overshoot   and  steady  state. 

3.  Effect    on  State   Variable   Due    to    Variatioi    in    C 

Figs. 3. 40  and  4.6  (^7%  and  Ai23)  show  that  <5rc  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected    in   the   steady  state. 

Figs. 3. 42  and  4.7  (Ag^and  A135)  show  that  &r,  is  not 
affected  in  the  rise  time  and  steady  state,  and  strongly 
affected  in  the    overshoot. 

Figs. 3.  46  and  4.8  (  A\o&  and  Aj^3  )  show  that  6y0  is 
little  affected  in  the  rise  time,  little  affected  in  the 
overshoot   and  not  affected    in   the  steady  state. 

Figs. 3. 48  and  4.9  (Aug  and  A153  )  show  that  by  is 
little  effected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected    in   the   steady  state. 

Figs. 3. 32  and  4.10  (^gand^i^  )  show  that  (^  is 
little  affected  in  the  rise  time,  stronglu  affected  in  the 
overshoot   and  not  affected    in   the   steady  state. 

Figs. 3. 28  and  4.11  (^,0  and  ^173 )  show  that  y  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot   and  not  affected    in   the   steady  state. 
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Figs. 3. 30  and  4.12  ( A2£  and  Al&3  )  show  that  b  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot   and  not  affected    in   the   steady  state. 

Figs.  3.  50  and  4.13  (  A (23  and  A  1^3)  show  that  £  is 
not  affected  in  the  rise  time  and  steady  state  and  strongly 
affected  in  the    overshoot. 

4.      Effect    on  State   Variable   Due    to    Variation    in   C  : 

Figs. 3. 40  and  4.6  (A74-  and  ^I2«q)  show  that  Src  is 
not  affected  in  the  rise  time  and  steady  state  and  little 
affected   in  the    overshoot. 

Figs. 3. 42  and  4.7  (Ae5  and  A134  *  show  that  Sr.  is 
not  affected  in  the  rise  time  and  steady  state  and  strongly 
affected   in  the   overshoot. 

Figs. 3. 46  and  4.8  (\i0$  and  A|^^  )  show  that  ofc  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot   and  not  affected    in  the   steady  state. 

Figs. 3. 48  and  4.9  (A(|$  and  \\S^)  sil°w  "that  Sy  is 
not  affected  in  the  rise  time  and  steady  state  and  strongly 
affected   in  the   overshoot. 

Figs. 3. 32  and  4.10  (  A3^and  ,X  K^  )  show  that  p>  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot   and  not  affected    in   the   steady  state. 

Figs. 3. 28  and  4.11  (Aisand  Ai7<a)  show  ttiat  r  is  not 
affected  in  the  rise  time  and  steady  state  and  strongly 
affected  in  the    overshoot. 

Figs. 3. 30  and  4.12  (A25  and  Al©4.  )  show  that  p  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot   and  not  affected    in  the   steady  state. 

Figs. 3. 50  and  4.13  (A12.5  and  A|g4)  show  that  <f>  is 
not  affected  in  the  rise  time  and  steady  state  and  strongly 
affected   in  the   overshoot. 
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5«      Effect    on  State   Variable   Dus    to    Variatioi    in   C : 

Figs. 3. 39  and  4.6  (^73and  A|2.6~)  show  that  Sfc^is 
little  affected  in  the  rise  time,  strongly  affscted  in  the 
overshoot   and  not  affected    in   the   steady  state. 

Figs. 3. 41  and  4.7  (Ap3and  A135  )  show  that  6&  is 
little  affected  in  the  risr  time,  stronglu  affscted  in  the 
overshoot   and  not   affected    in   the   steady  state. 

Figs. 3. 45  and  4.8  (A^o^  and  ^l-^S  )  show  that  6yc  is 
little  affected  in  the  rise  time,  strongly  affscted  in  the 
overshoot   and  not  affected    in   the   steady  state. 

Figs. 3. 47  and  4.9  (  An  3  and  A^j)  show  that  £>y  is 
little  affected  in  the  rise  time,  strongly  affscted  in  the 
overshoot  and  not  affected    in  the   steady   state. 

Figs. 3. 31  and  4.10  O33  and  Ai6S  )  sho*  that  fi  is 
little  affected  in  the  rise  time,  stongly  affected  in  the 
overshoot   and  not  affected    in   the   steady  state. 

Figs. 3. 27  and  4.11  (A13  and  A  17S  )  show  that  V  is 
little  affected  in  the  rise  time,  strongly  affscted  in  the 
overshoot   and  not  affected    in   the  steady  state. 

Figs. 3. 29  and  4.12  (/\2^and  A|g5  )  show  that  \>  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot   and  not  affected    in   the   steady  state. 

Figs.  3.  49  and  4. 13  ( A  123  and  Ai9£)  show  that  £  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot   and  not  affected    in   the   steady   state. 

C.       EFFECT    OF   TRAJECTORY   SENSITIVITY    FUNCTIONS 

1 .      Oncoupled  Pitch    Autopilot 

The  trajectory  sensitivity  functions  ^11,  /112, 
>L13,  and  ,114  in  the  linear  case  and  the  correspondent 
trajectory  sensitivity  functions  151  and  152  ia  the  nonli- 
near case  show  strong  effect  on  ths  rise  time  aad  overshoot 
of   the  state  variable  q. 
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The  trajectory  sensitivity  functions  /L2  3  and  ,124 
in  the  linear  case  and  the  correspondent  trajectory  sensi- 
tivity function  16  1  in  the  nonlinear  case  show  strong 
effect  on  the  rise   time   of    the   stata    variable    c*    . 

The  trajectory  sensitivity  functions  ,13  3  and  ,134 
in  the  linear  case  and  the  correspondent  trajectory  sensi- 
tivity function  /14  1  in  the  nonlinear  case  show  strong 
effect  on  the  overshoot  and  steady  state  of  the  state  vari- 
able  <5p  . 

The  trajectory  sensitivity  functions  ,15  3  and  /164 
in  the  linear  case  and  the  correspondent  sensitivity  trajec- 
tory function  .131  in  the  nonlinear  case  show  strong  effect 
on   the  overshoot   and   steady    state   of    the   state  variable  $p  . 

2.      Coupled    Roll -Yaw   Autopilot 

The  trajectory  sensitivity  functions  ,113, /115,  and 
,118  in  the  linear  case  and  the  correspondent  trajectory 
sensitivity  functions  1.175,  /1174,  and/1173  in  the  nonlinear 
case  show  strong  effect  on  the  overshoot  of  the  state  vari- 
able   r. 

The  trajectory  sensitivity  functions  ,123  ,  /i25,  and 
128  in  the  linear  case  and  the  correspondent  trajectory 
sensitivity  functions  ,1185,  1.184,  and  ^L 1 83  in  the  nonlinear 
case  show  strong  effect  on  the  overshoot  of  the  state  vari- 
able   p. 

The  trajectory  sensitivity  functions  ,13 3  ,  135,  and 
138  in  the  linear  case  and  the  correspondent  trajectory 
sensitivity  functions  1.165,  ,1164,  and/1163  in  tie  nonlinear 
case  show  strong  effect  on  the  overshoot  of  the  state  vari- 
able   fh    . 

The  trajectory  sensitivity  functions  /173  ,  .175,  and 
178  in  the  linear  case  and  the  correspondent  trajectory 
sensitivity  functions  1.125,  /1124,  and /1 123  in  tae  nonlinear 
case  show  little  effect  on  the  overshoot  of  the  state  vari- 
able   So   . 

166 


The  trajectory  sensitivity  functions  ,183,  /185,  and 
,188  in  the  linear  case  and  the  correspondent  trajectory 
sensitivity  functions  1135,  ,1134,  ind/1133  in  tae  nonlinear 
case  show  strong  effect  on  the  overshoot  of  the  state  vari- 
able   5r  - 

The  trajectory  sensitivity  functions  /i1  03  and  ^.105 
in  the  linear  case  and  tha  correspondent  trajectory  sensi- 
tivity functions  ,1145  and  ,1144  in  the  nonlinear  case  show 
strong  effect  on  the  overshoot  of  the  state  variable  oYc  • 
The  trajectory  sensitivity  function  ,11 08  in  the  linear  case 
and  the  correspondent  trajectory  sensitivity  f  met  ion  /1143 
in  the  nonlinear  case  show  little  effect  on  the  overshoot  of 
the   state  variable  Sy0  • 

The  trajectory  sensitivity  functions  ,1113,  ,11 15, 
and  ,1118  in  the  linear  case  and  the  correspondent  trajec- 
tory sensitivity  functions  ,11 55,  ,1154,  and  ,1153  in  the 
nonlinear  case  show  strong  effect  on  the  overshoot  of  the 
state   variable    Sy  . 

The  trajectory  sensitivity  functions  /1123,  ,1125, 
and  1128  in  the  linear  case  and  the  correspondent  trajec- 
tory sensitivity  functions  ,1195,  ,1194,  and  A.  193  in  the 
nonlinear  case  show  strong  effect  on  the  overshoot  of  the 
state  variable    <£     . 

The  trajectory  sensitivity  functions  have  the 
following  range    of    values   : 

Minimum  Maximum 

Linear  -  385.25  632.35 

Nonlinear  -  642.9  1772.9 


D.       GENERAL   CONCLUSIONS 

The      parameter   C_         strongly      affect      the  o/ershoot      in 
almost  all    case. 
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The    parameter     C..        little   affect      the   overshoot      in  all 

Mot- 
case. 

The     parameter    C0  strongly      affect      the   overshoot     in 

almost  all   case. 

The      parameter    C„  strongly      affect      the   overshoot     in 

almost   all   case. 

The      parameter    C^  strongly      affect      the  overshoot     in 

almost  all   case. 


E.       RECOMMENDATIONS 

It  was  not  possible  to  run  a  CSMP  computer  program  using 
all  sta-ce  variable  (19)  and  parameters  (10)  in  the  nonlinear 
case  due  to  the  work  area  available.  The  above  restriction 
occurred  just  when  the  actual  and  nominal  output  were 
required  to  be  printed.  The  computational  time  in  this  case 
may  be  reduced  using  Fortran  subroutines  imbedded  in  the 
CSMP  program  in  order  to  calculate  the  necessary  parameter 
derivatives. 

In  the  linear  cases  the  computational  time  can  be 
reduced  by  means  of  using  the  method  of  Sensitivity  Points 
that  uses  one  model  for  all  parameters  instead  of  using  one 
model    for  each    parameter. 

The  present  analysis  can  be  repeated  for  the  case  of  the 
Circular2  Airframe  given  in  the  [Bef.  2  ],  and  comparisons 
between   both  airframes  can    be    performed. 

Future  study  can  be  made  modeling  the  system  in  the 
frequency  domain   where  the    "root   sensitivity"    is    analysed. 

Influence  of  parameter  variations  on  miss  distance  can 
be  analysed  using  an  augumented  system  and  a  scenario  of 
reference. 


2The      current      CSMP      programs        in      the      appendices      are 
prepared  to   run    the    Circular  case   just   by   inserting   the  data 


of  the   correspondent   airframe. 
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APPENDIX    A 
MISSILE   SIZING,     MASS    PROPERTIES    AND    AERODYNAMIC    DATA 

A.  INTRODUCTION 

In     this  appendix,  one     considers      a   model      that      was 

assumed  to  be  1/6  scale  of  the  actual  elliptical  missile 
configuration  as  given  in  Fig.  A.  1. 

Missile  canf iguration  was  sized  to  provide  realistc  mass 
properties    needed  for  this    study. 

Since  the  main  purpose  of  this  work  was  to  perform  the 
parameters  sensitivity  analysis  applied  to  a  3ank-to-Turn 
missile,  no  effort  was  expended  oa  a  detailed  design  and 
analysis  of  the  various  parts  envolved,  which  is  given  in 
[Ref.  2].  Some  of  figures  given  in  this  appendix  were  taken 
from  the  [Ref,  2]  fcr  easy  visualization  of  the  system  in 
study. 

B.  GEOMETRIC  AND  MASS  PROPERTIES  OF  MISSILE  CONFIGURATION 

One  can  find  in  Table  VI  the  properties  used  in  the 
development  of  the  equations  for  applying  the  Parameters 
Sensitivity   Analysis. 

Fig. A. 2  shows  the  aerodynamic  sign  convention,  nomencla- 
ture and  body-fixed  axes  where  the  following  three  assump- 
tions  were    made: 

-  Plane  x    -  zi    of    Fig. A. 2   is   the   maneuver   plane. 

-  Missile    is   trimmed    in   pitch    (i.e.    ,    My=0    ,    at   fixed 

values   of   q,<x,    and   £p  )  . 
-Missile    roll    rate   is   constant. 

Fig.    A.  3  shows    a    block    diagram   of   a    BTT   autopilot. 
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Figure  A.I   Hodel  of  elliptical  configuration  (  1/6  scale  ) 
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TABLE   71 

Geometric    and   Mass  Properties 
Configuration. 


of   Missile 


WEIGHT  (lbs) 24  75 

IXx(slug   sq    ft) 110 

Iy/(slug   sq    ft) 790 

(slug   sq    ft) 953 


zz. 


LENGTH    (in) 


1  68 


CENTER   OF   GRAVITY 100.8(0.6    1) 

distance   from   nose    (in) 


MAX.     MAJOR    AXIS     (in) 
MAX.     MINOR    AXIS     (in) 


41.7 
13.86 


Inertial  acceleration  commands  are  applied  in  pDlar  coordi- 
nates (i.e.,  magnitude  of  command  ( T|c  )  applied  to  the 
pitch  autopilot  and  the  direction  (  <fcc )  is  applied  to  the 
roll    autopilot)  . 

The  yaw  autopilot  is  slaved  to  the  roll  autopilot  to  mini- 
mize sideslip  angle  by  coordinating  the  missile  yaw  and  roll 
motion. 

Achieved  maneuver  plane  or  inertial  acceleration  in  rectan- 
gular coordinates  (i.e.,  *L  and  "*]yv/J  ^s  determined  by 
resolving  achieved  body-fixsd  accelerations 
(i.e.,1l  and  >|  )  through  missile  roll  angle(<|>  )  (i.e., 
Euler  angles  ©■  and  <f>  are  assumed  to  be  sufficiently 
small) . 


171 


Figure  A.  2    Aerodynamic  Sign  Convention 
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APPENDIX    B 
LINEAR    SYSTEM    H3DELS 

A.  INTRODUCTION 

The  following  appendix  addresses  the  informations  of  the 
linear  uncoupled  pitch  channel  autopilot  and  the  roll-yaw 
coupled  channel  autopilot  of  a  elliptical  airframe  taken 
from    [Ref.    2]. 

The  linear  time  domain  analysis  of  the  CBTT  autopilot  was 
assumed  that  the  missile  is  initially  in  the  desired 
maneuver  plane  and  trimmed  at  ten  degrees  anglas-of-attack 
(i.e.,       the  squilibrium  or      trim   angle-of-attack  in  the 

model  of  Fig. A.  1  equals  10  degrees  and  the  equilibrium  roll 
rate  Pe  in  the  models  of  Fig. A. 3  sgual  zero). When  Pe  =  0, 
Qe  (i.e.,  equilibrium  pitch  rate)  has  been  formed  to  have 
negligible  influence  in  the  lateral  model  compared 
to    o<      and   was    therefore   set   equal  to  zero. 

B.  UNCOUPLED  PITCH  CHANNEL  AUTOPILOT 

A  general  block  diagram  of  an  uncoupled  pitch  channel 
autopilot  is  shown    in  Fig.    B.1  . 

A  normal  acceleration  command  (V)z,  g's)  is  applied  to 
the  pitch  control  law  which  uses  measurements  of  missile 
body  pitch  angular  rate  (q)  and  pitch  normal  acceleration 
(*l2)  to  determine  the  required  actuator  command  (Sp  )  .  The 
actuator  is  modeled  as  a  first-order  lag  at  30Hz  (188.4 
rad/sec)  . 
The  dynamic  model  is  linearized  about  a  trim  angle-of- 
attack  as  described    in  appendix   A. 

From  the  block  diagram  of  Fig. B.1,  the  parameters  of 
interest      for   the     present    analysis      are   given     respectively 


as 


^±i       »v#ii„»Lrtt_#V» 


N 


Ot 


"V     "*<*    '    ^«p 


174 


f  I 

^_ 

■*r— 

■ 

"*~ 

• 

I 

0 

*/> 

JC 

OQ 

*"* 

Figure  B. 1 


Uncoupled  Pitch  Channel  Autopilot, 
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For  ease  notation  the  correspondence  of  state  variables  and 
constants  are  given  in  Table  VII 

Performing  block  diagram  manipulations  see  [Bef.  5],  one 
yields  to  the  following  state  variable  equations: 

Aerodynamic  equations 
X1  =  C2  A3  X2  ♦  C2  A4  X3  (B.1) 


• 


X2  =  X1  -  K  C1  A1  X2  -  K  C1  A2  X3  (B.2) 

Control    law   equations 
X3    =  -  C3   X3    ♦   C3    Conv   X6  (B.3) 


14    s   -   CI   CH    A1    X2    -   C\    C4    X3    -    C4    XH  (B.4) 


• 


X5    =   C7    XU   -    C5   X5    -   C6    NZC  (B.5) 

Actuator    equation 
X6    =   C    E    A7    X3    -    6    X6    ♦    C    E    A8    X8    ♦    C    E    A6    X11  (B.6) 


C.       COUPLED   ROLL-JAW    CHANNEL   AUTOPILOT 

A  general     block   diagram  of     a  coupled     roll- yaw  channel 
autopilot  is  shown    in  Fig. B.2. 
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Figure  B.2        Coupled    Boll- Yaw  Channel  Autopilot. 
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The   block   diagram    (Fig.B.2)       shows,      eight  parameters   of 
interest  that   are: 

A1    =    CYSy'    A2   =    CYp     '    A3    =   C*1Sr'    AU    =    Cn^     ' 

For  ease  of     notation  the   correspondence   of     state   vari- 
ables and  constants    is  given  in  Table   VIII. 

Eqns.B.7  through   B.17  give      the   state   variable   equations 
of   the   coupled    roll-yaw  channel   autopilot: 
Aerodynamic   equations 

X1    =   C  Conv(AU    X3    ♦   A3    X8    ♦    A5    X11)  (B.7) 


X2    =   C   Conv(A7    X3    ♦    A8    X8    +    A6    X11)  (B.8) 


X3    =  -   X1    -    (alpha/Conv)     X2  (B.9) 

♦  KB    A    A2  X3    ♦    KB    A   A1    X11 

X4    =   -8    XU  -17.6    X12   ♦    17.6    PHC  (B.  10) 

Control  law  equations    (roll) 
X5    =   -     (0.755/Conv)    X2   -    C   D    A7    X3  (B.  11) 

♦  (0.755   -    8    D)X4    -   5    X5    -    C    D    A3    X8    -    C    D    A6    X11 
-     17.6    D    X12    +17.6   PHC 
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X6    =   C    E    A7    X3   -    6    X6    +    C    E    A8    X8    ♦    C    E    A6    X11  (B.12) 


X7    -   -    F    KC    (0.755/Conv)     X2    -    (D   +    E)     F    KC    A7    X  3        (B.  13) 

♦  F    KC(0.755    -   8    D)    M   ♦    KC(15  -    5    F)     X5 
«-KC(6F   -    15)     X6    -    15    X7    -    (D    +  E)     F    KC    C    A8    X3 
-     (D    +    E)    F    KC  C    46   X11     -    F    KC  D    17.6    X12 

♦  F    KC    D    17.6    PHC 

X8    =    188.4  Conv   X7    -    188.4    X8  (B.14) 

Control  law   equations    (yaw) 
X9    ■    (K1    A   A2/T1)     X3    -   X9/T1    +     (K1     A    a1/Tl)     X11  (B.15) 


X10   =     (K2/Conv)    X1    -    (K2    H   ALPHAB)  /Conv    X2  (B.  16) 

♦  (K2/10)((K1    A    A2)/T1    ♦    B    A4    -    a     ALPHAB    C    A7    X3 

♦  (K2/10)  (B    A3   -    H    ALPHAB   C    A8    X3    ♦    K2(1    -     (1/10    T1)     X9 

♦  (K2/10) ((K1    A    A1)/T1)     ♦    B    A5    -    H    ALPHAB    C    A6)     X11 

Actuator  equations 

X11    =    188.4   Conv    X10   -    188.4    X11  (B.  17) 


X12   =    X2/Coav  (B.  18) 
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D.       AERODYHAHIC    DATA-LINEAR    APPROXIMATION 

A  linear  approach  was  used  in  the  design  aad  stability 
analysis  of  the  aultopilots  of  tha  pitch,  yaw,  and  roll 
channels,  both  uncoupled  and  coupled.  According,  a  linear 
approximation  of  the  aerodynamic  dsrivatives  at  m=3.95  was 
provided,  aboat  which  the  system  could  be  perturbed.  These 
linearized  aerodynamic  derivatives  are  given  in  [  Ref .  2]  and 
are    presented   here    in  Table    IX. 


180 


TABLE   ¥11 
Correspondence  of   Symbols    (Uncoupled   Pitch   Autopilot) 


X 
Y 

cm 


sP 


qs/d 

57,3qsd/I 

188.4 

150.0 

6.0 

0.53544 

0.48 

0.38375 

3.07 

0.48 

57.3 


XI 

X2 
x3 
X4 
x5 

X6 
A1 
A2 
A3 
A4 
C1 
C2 
C3 
C4 
C5 
C6 
CI 
C8 
C9 
K 
Conv 
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TABLE   7III 

Correspondence  of   Symbols    (Couplsd   Roll-Yaw  Autopilot) 


r  X1 

p  X2 

£  X3 

X  X4 

Y1  X5 

XI  X6 

£ftc  x7 

6r  *8 

y  X9 

&Yc  X1° 

Sy  X11 

$  X12 

CV/.  A1 

CY^  k2 

Cn^  A3 

Cft  AU 

O  A5 

%K  A6 

C^  A? 

%* &8 


3. 1416  S 

57.3  Conv 

qs/d  A 

qsd1/I22  B 

qsd1/IXx  C 

0.05033  D 

0.078  E 

0.25  F 

*  €  ALPHA 

°<>e  /57.3  ALPHAB 

0.48  KB 

0.25  T1 

0.84  K1 

6.08  K2 

1*0  KYP 

4.17  KC 

0,48  K 

2.0  D1 

q  QB 


182 


TABLE   IX 
Linearized   Aerodynamic   Derivatives 


Hc< 


N 


S? 


*1o< 


m 


*p 


•x  =   100 

-  0.05a 

♦  0. 02a 

-  0. 027 

♦  0. 015 

-  0. 039 

-  0.  010 

-  0. 006 

♦  0. 014 

♦  0.  023 

♦  0.  220 
+    0.  020 

♦  0. 0137 

-  0. 055 
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APPEND*!   S 
NONLINEAR   SYSTEM    MODEL 

A.       INTRODUCTION 

A  general  block  diagram  of  the  nonlinear  moiel  is  shown 
in  Fig. C.1  The  three  dimensional  nonlinear  aerodynamic  model 
presented  hers  is  for  the  same  conditions  used  ia  the  linear 
model   shown  in    appendix   B.  The   same   flight   condition  used 

for  the  linear  case  is  the  same  for  the  nonlinear  model  (i. 
e.  ,  60  kft  altitude,  Mach  3.95).  The  control  Laws  are  the 
same  used  for  the  linear  models  except  for  a  minor  modifica- 
tion to  the  coordinating  branch  dependence  on  angle-of- 
attack  and  also  the  inclusion  of  anti-gravity  bias  as  stated 
in  [Ref.    2]. 


B.       CONTROL   LAWS 

The  control  laws  used  for  the  following  nonlinear  3-D 
studies  were  the  same  as  those  used  for  the  linear  models, 
except  for  the  gain  o<  shown  in  the  bold  line  of  coordina- 
tion branch  in  Fig.  C.1.  The  new  gain  oc  is  held  constant 
at  one  degree  magnitude  for  angle-of-attack  less  than  one 
degree     positive.  For  angle-of-attack      greater      than     one 

degree  positive,  the  gain  °<  is  equal  to  the  angle-of- 
attack.  This  keeps  coordination  for  very  small  angle-of- 
attack. 

Gravity  effects  were  not  included  in  the  linear  models 
in  appendix  B,  because  it  was  assumed  to  have  a  negligible 
influence  on  autopilot  stability  and  response  for  pertuba- 
tions   about     a    missile     trim  condition.  However,      gravity 
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Figure  C. 1    Nonlinear  Model. 
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effects  were  included  in  the  following  nonlinear  model  where 
the  missile  body-fixed  yaw  axis  will  be  subjected  to  the 
full  force  of  gravity  and  may,  therefore,  have  a  significant 
influence  on  sideslip. 

The  gravity  effects  are  compensated  for  by  pitch  and 
roll  motion  of  the  missile  which  have  less  influence  on 
sideslip  than  yaw   motion.  In  inertial   rectangilar   coordi- 

nates one  has: 

^|  c  =acceleration   ccmmand   in   inertial    zy  direction 


'V  =  ^7     "    Cos<e>  <c-1> 

where 
^2.       =  guidance  command(    gfs   ) 
-Cos  e   =  anti-gravity  bias    command    (    g*s   ) 
#v\vc    =   acceleration  command   in   the    inertial    y      direction 
guidance   command  (  g's  ) 
There   is  no    gravity   effect   in   Y.    direction.        In   order  to 
assure   an     anti-gravity   bias     of    just   one     g,      oa e     needs   to 
modify  the   anti-gravity  command  as  follows: 

elliptical    airframe   inertial   acceleration  command 
*\      -  Nzc-0.913cos(<(>) 


C.       PEBFOHANCE    COMMANDS 

The   commands        2      gees      is   first      applied   in      the   0°     or 
upward  direction     at      2     seconds.  Since    both     the   missile 

roll  angle  and  roll  angle  command  are  at  zero  degrees,  there 
is  no  roll  motion  and  the  missile  turns  upward  as  a 
skid-to-turn  controlled  missile.  At  5  seconds,  a  second 
2     g's   inertial    guidance  command      is    applied    in   the   downward 
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or  180°  direction.  The  missile  is  commanded  to  roll  through 
180°  while  moving  in  a  coordinated  lanner  in  yaw  and  roll  to 
minimize  sideslip  angle  and  prevent  or  minimize  negative 
angle -of -at tack. 

The  response  of  achieved  maneuver  and  cross-plane  accel- 
erations during  the  first  2  seconds  are  due  to  initial 
conditions,      gravity     and  anti-gravity     bias   effects.  The 

initial  conditions  were  added  to  minimize  the  transients 
which  result  when  gravity  bias  commands  the  autopilot  for 
constant  alritude  missile  flight.  The  initial  conditions 
were : 

<*      ■   angle-of-attack   =  2.41   deg 
dp     =   pitch  tail    angle   ■    0.658    deg 
©      =    pitch  Euler    angle  =    3.65    deg 

output    of  pitch   acceleration    feedback   lag    =   -1.0   g*s 
pitch   control  law  acceleration   error 
lag  prior  to   dc   gain  =  -    0.3105 
6|>     =   pitch   actuator  command   =   0.658   deg 

The  achieved  maneuver  plane  accsleration  (Egn.C.2)  is 
calculated    from    the    body-fixed  acceleration    T)       ,    V\  ,    and 

the   roll  angle    Q  as    follows: 


*U=    lie005^    *  *W6    Svn  * 


(C.2) 


During  the  first  command,  achisvsd  body-fixei  yaw  accel- 
eration (l^yo)  and  missile  roll  angle  are  egual  to  zero, 
because     the  roll     channel      is      not   commanded.  Therefore, 

achieved  maneuver  plane  acceleration  is  egual  to  the  body- 
fixed  acceleration  V[  z  .  During  the  second  command  the 
missile  roll  angle  has  the  same  roLl  angle  response  of  the 
roll    autopilot. 
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D.       AERODYNAMIC    DATA    -   NONLINEAR    REPRESENTATION 

The  aerodynamic  data  are  given  in  [Ref.  2].  The  entire 
study  was  conducted  using  M=3.95  and  the  aerodynamic  coeffi- 
cients  are   based   on    a  body- fixed   axis   system   of   Fig.B.2. 

For  reference,  the  normal  force  and  pitching  moment 
plots  given  in  appendix  B  of  [Ref.  2  ]  at  M=3.  95  were  repro- 
duced in  Table  X  and  XI,  which  give  just  the  values  for 
points  of  interest.  These  tables  present,  too,  the  corre- 
spondent derivatives  that  were  found  by  simulation  using  a 
fortran  program  shown  in  appendix  F.  The  aerodyiamic  deriv- 
ative of  Cjj  ,  C-  ,  and  Cy  with  respect  to  sideslip 
angle   (b     ,        yaw     control  6y  ,  and     roll      control  S^     are 

presented  in  Table  XII,  XIII  and  XIV  which  give  just  values 
for  points  of  interest  and  includes  the  correspoadent  deriv- 
atives that  were  found  by  simulatioa  using  a  fortran  program 
given   in   appendix  F. 

The  block  diagram  of  Fig.C.1  shows  10  parameter  of 
interest  that  are   : 

A1    =    Cm(*c,<fp)     ,     A2   =   C^(k  ,fy)     .    A3    =    C$       (*)     , 

A4    =    Cnc    (oc)        ,     A5   =   C-        («.)       ,    A6    =    CYc    (  <*. )     , 

A7    =    Cfi        <o<)       ,     A8   =   C«       (X)     ,       A9    =   C         {<*  )     , 

and  AA  =  CY.    (<*)  • 

Using  the  same  procedure  as  given  in  appendix  B  one  can 
have  the  Tables  XV  and  XVI  that  give  the  correspondence  of 
symbols  for  the  nonlinear  model.  Egns.C.3  through  C.21  give 
the   state  variable    equations  of   the   nonlinear  moiel. 
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X1  ■  -  C»  X1  -  C02  CU  CN(X6,X4)  (C.3) 


X2  s  C7  II  -  C5  X2  ♦  C6  C19  Cos(X17)  -  C6  NZC       (C.4) 


X3  =  -  C7  C8  X1  ♦  (C5  C8  -  C9)  X2  ♦  (C9/Conv)  K  5    (C.5) 

♦    C6    C8    NZC 
XU    =   -    C3    X4    ♦    C3    Conv    X3  (C.6) 


4 

X5    =    (X17   X18)/Conv  ♦    C1    Conv    CM(X6,X4)  (C.7) 


X6    =    X5    -    KB    CM(X6,X4)     -     (X16    X18>/Conv  (C.8) 


X7    =   X5   Cos  (X  19) /Conv    -   X17    Sin  (X1  9) /Conv  (C.9) 


X8    =    -    C10    X8    +    (C12   -    C11    C13)     X10  (C.  10) 

♦  C11    C1U    PHC    -   C11   C1U    X19    -    CO    C1l(c„     (X6)     X1  5 

♦  C0„     (x6)     X13)-    C12    X18/Conv  T 

"or, 

X9    =   -    X9/T1    ♦     (K1/T1)     C02(C       (X6)     X16  (C.11) 

Y(* 

♦  C       (X6)    X15) 


189 


X10    =    -    C13    X10    ♦    C14    PHC    -    C14    X19  (C.  12) 


X11    =    -    C15    X11    «■    C16    CO(C?      (x6)     X16  (C.  13) 

+    C       (X6)    X15    ♦    C       (X6)     X13) 
X12    =    -    C7    X12    ♦    K  (C17    -    C10/C18)     X18  (C.  14) 

♦  (K/C18)  (312    -    C11    C13)      X10    ♦    K    C11     (C1tt/C18|      PHC 

-  K    C11    C14/C18    X19    -    (K/C18)     CO    C        (X6)(C11 

♦  C16)    X16    -     (K/C18)     CO    CJ(3(X6)(C11    ♦    C16)     X15 

-  (K/C18)    30    C^(X6)     (C1  1    ♦    C16)     X13 

-  (K    C12/CDI1V)     X18    ♦     (C15/C18    -    C17)     XI  1 

X13    =   -   C3  X13   ♦    C3  Conv    X12  (C.  15) 


X14    =     (K2/10)   (-    X9/T1    ♦    K1    T1    C02     (C       (X6)     X16  (C. 16) 

♦  CXs  (X6)    X15)    -     (X15    X1  8)/Conv2(X5    -    KB    C02  CN(X6,XU) 

-  (X16   X18)/Conv    X18   -    (KYP/Conv)  (X5    -    KB    C02    CN(X6,X4) 

-  (X16   X18)/Conv    X18    -    (KYP/Conv)     CO    XX6     (C^    (X6)     X16 

-  Cjj    (X6)    X15    ♦    CSs   (X6)     X13)     ♦    (K2/Conv)     X17 

♦  K2    X9    -    (K2    KYP/Conv2)     XX6    X18 


X15    =    -   C3   X15    ♦    C3   Conv    X14  (C. 17) 


X16    =    KB    C02  (C       (X6)     X16     ♦    C       (X6)     X15  (C.18) 

T(>  iff 

♦     (X16    X18)  /Conv    -  X17 
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X17    =   -     (X5   X18)/Ccnv   ♦   C01    Conv(3n     (X6)     X16 


(C.  19) 


♦    C^(X6)    X13) 


% 


X18   =    CO    Con 


v  (C«     (X6)     X16    +    C0 


(X6>     X15 


(C.  20) 


♦    C       (X6)    X13) 


X19    =    X18/CDnv 


(C.  21) 
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TABLE   X 
Normal  Porce   Coefficient  and  Derivatives  M=3.95 


cNc*,V> 


"^C^U^p)         ICm^p) 


"~bu 


!>6 


-4.00C0 

-10.000C 

-1.C500 

0.1842 

0.0308 

-4. COCO 

0.0000 

-0.8000 

0.1453 

0.0208 

-4.C0CC 

5.000C 

-0.7000 

0.2496 

0.0196 

-4.0000 

10.0000 

-0.6000 

0.2141 

0.0208 

0.00CG 

-10.0000 

-0.2500 

0 .2110 

-0.0332 

0.00C0 

0.0000 

-o.caoo 

0.20  74 

0.044  8 

0.00  00 

5.0000 

0.1300 

0.1827 

0.0336 

0.0000 

10.0O0C 

0.2J00 

0.1929 

-0.0112 

4.0000 

-10.0000 

0.6000 

0.2092 

0.020C 

4.00CC 

0.0000 

0.8000 

3.2253 

0.0200 

4.C3CC 

5.O00C 

C.9000 

0.2196 

0.0200 

4.0000 

10.000C 

l.COOO 

0.2141 

0.0200 

8. COCO 

-10.000C 

1.4500 

0.2272 

0.0042 

8.0000 

0.000C 

1.7C00 

0.2265 

0.0342 

8.0000 

5.0000 

1.8500 

0.2288 

0.0229 

8.00  00 

10.000C 

1.9000 

0.2256 

-0.0058 

12.0000 

-10.0O0C 

2.4000 

0.2321 

0.0200 

12.0000 

0.0000 

2.6000 

0.2188 

0.0200 

12. COCO 

5.000C 

2.7000 

0.2153 

0.0200 

12.0000 

10.OOOC 

2.8000 

0.2335 

0.0200 

16.C0CC 

-10.000C 

3.3000 

0 .2319 

0.0083 

16.0000 

0.000C 

3.50C0 

0.2482 

0.0283 

16.0000 

5.0000 

3.6500 

0.2599 

0.0308 

16.00CC 

10.000C 

3.6000 

0.2654 

0.0283 

20.0000 

-10.0000 

4.3000 

0.2653 

0.0300 

20. COCO 

o.oooc 

4.6000 

0.2884 

0.0300 

20.0000 

5.0000 

4.7500 

0.2825 

0.030C 

20.0000 

10.0000 

4.9000 

0  .2798 

0.0300 

24.0000 

-10.000C 

5.4GC0 

0.2819 

0.0300 

24.0000 

o.oooc 

5.7000 

0.2482 

0.0300 

24.C0CC 

5.000C 

5. £500 

0  .2599 

0.0300 

24.0000 

10.000C 

6.C0OO 

0.2654 

0.0300 
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TABLE   XI 
Pitching  Homent   Coefficient   and    Derivatives    M=3.95 


tnW 


"~<)Cm(°Up)    ^Cmi^p) 


D* 


TW, 


-4.0O0C 

-10.000C 

0.4850 

0.0130 

-0.0159 

-4.0000 

-5.0000 

0.2800 

-0.0200 

-0.C61C 

-4.0000 

0.0000 

-0.C750 

0  .0190 

-0.0759 

-4.00C0 

10.000C 

-Q.£300 

0.C197 

-0.0149 

0.0000 

-10.0000 

0.5500 

0.0185 

-0.0550 

o.oocc 

-5.000C 

0.2750 

0.0131 

-0.0550 

0.0000 

0.0000 

-o.cooo 

0.0186 

-0.055C 

0.0000 

10.0000 

-0.5500 

0.0195 

-0.0550 

4.00CC 

-10.000C 

0.6250 

0.0180 

-0.0551 

4.0000 

-5.0000 

0.2500 

0.0200 

-0.055C 

4.0000 

0.0000 

0.C750 

0.0190 

-0.0551 

4.0000 

10.000C 

-0.4800 

0.0147 

-0.0561 

8.0000 

-10.0000 

0.7100 

0.02  94 

-0.0562 

8.00C0 

-5.000C 

0.4300 

0.0232 

-0.0558 

8.0000 

o.oooc 

0.1500 

0.0178 

-0.0563 

8.0000 

10.0000 

-0.4300 

0.0118 

-0.0603 

12.000C 

-10.000C 

0.6500 

0.0331 

-0.0640 

12.0000 

-5.0000 

0.5250 

0.0184 

-0.0655 

12.0000 

o.oooc 

0.2000 

0  .0034 

-0.064Q 

12.0000 

10.0000 

-0.2900 

0.0057 

-0.0520 

16.C000 

-10.0000 

0.9500 

0.0182 

-0.0725 

16.0000 

-5.000C 

0.5750 

0.0120 

-0.0763 

16.0000 

0.0000 

0.2000 

0.0059 

-0.0725 

16.00CC 

10.000C 

-0.4000 

-0.0122 

-0.0425 

20.0300 

-10.000C 

1.CC0O 

0.0065 

-0.0756 

20.0000 

-5.0000 

0.6220 

0.0065 

-0.0755 

20.03CC 

o.oooc 

0.2450 

0.0067 

-0.0752 

20.0000 

10.0000 

-0.5000 

-0.0395 

-0.0736 

24.00CC 

-10.0000 

1.C000 

-0  .006  8 

-0.0841 

24.0000 

-5.000C 

0.5870 

-0.0290 

-0.0818 

24.0000 

0.0000 

0.1750 

-0.0516 

-0.0837 

24.00CC 

10.000C 

-0.7300 

-0.0772 

-0.1001 
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1 — 

TABLE   III 

Sideslip 

Derivatives  ani 

Derivatives 

<* 

-4,0000 

0.0110 

-0.0027 

o.oooc 

0.0000 

-0.GO27 

4.00CC 

-0. C110 

-C.0022 

5.20C0 

-0.0135 

-0.0024 

8.C0CC 

-0.C220 

-0.0027 

10.0300 

-0. C270 

-C.0021 

12.6000 

-0.C310 

-0.0015 

16.00CC 

-0. C360 

-0.0  014 

20.0000 

-0.0410 

-0.0011 

<* 

v 

^(V*> 

~L)  c< 

-4.00C0 

0.C240 

CO  000 

0.0000 

0. C240 

C.  0  0  00 

4.0000 

0.C240 

G.00C0 

5.20CC 

0.0240 

0.0000 

8.00C0 

0. C240 

C.00C0 

10.0000 

0.C240 

C.0005 

12.60CC 

0.C270 

0.0009 

16.0000 

0. C290 

C.  0  0C7 

20. COCO 

0.0320 

0.0008 

<* 

Cy    U) 

-4.00CC 

-0. C490 

0.0  024 

0.00  00 

-0.0440 

0.OOC1 

4.CO0C 

-0.0480 

-0.0009 

5.20CC 

-0.  C490 

-C.0  0C9 

8.0000 

-0.C520 

-0.0013 

10.0000 

-0. C550 

-0.0013 

12.6000 

-0. C580 

-C.0010 

16.0000 

-0.C610 

-0.0009 

20.00C0 

-0.C650 

-0.0011 
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TABLE   XIII 

- 1 

Taw  Control 

Derivatives  atid 

Derivatives 

<x 

Co   <*) 



**x 

DoL 

-4.0000 

0.C025 

-C.  00C7 

0.0000 

0.0000 

-0.0006 

4.00CC 

-0. C020 

-0.0  014 

5.20CO 

-0.0040 

-C.0015 

8.00C0 

-0.C075 

-0.0013 

10. COCO 

-0.  C100 

-C.0014 

12.6000 

-0.0140 

-C.0015 

16. COCO 

-0. C190 

-0.0014 

20.0000 

-0.C240 

-0.0011 

o< 

-4. COCO 

-0.C400 

0.0  0  00 

0.0000 

-0.0400 

cooco 

4. OOOC 

-0.C400 

c.ooco 

5.20CC 

-0. C400 

0.0  000 

8.0000 

-0.C400 

cooco 

10.00CC 

-0.0400 

c.ooco 

12.60CC 

-0. C400 

-C.0  0C3 

16.0000 

-0.G420 

-0.00C7 

20.00CC 

-0. C453 

-0.0010 

OC 

&v    U.) 

^  C  Y^y^y 

Y*y 

^c* 

-4.00CC 

0. C150 

0.0  000 

0.0000 

0.0150 

cooco 

4. OOOC 

0.C150 

0.0000 

5.20CO 

0. C150 

cooco 

8. OOOC 

0.C150 

c.ooco 

10.C0CC 

0.C150 

0.0000 

12.60CC 

0.C150 

C0  0C3 

16.0000 

0.C170 

C.0006 

20.00CC 

0. C195 

0.0  006 
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TABLE    XIV 
Boll   Control  Derivatives  and  Derivatives 


^  dU^l 


<x> 

-7W 

-4.0000 

0.0232 

CO  000 

O.COCC 

0.0232 

0.0000 

4.0000 

0.C232 

cooco 

5. 2000 

0.C232 

c.  ooco 

8.0000 

0.0232 

0.0000 

10.0000 

0.  0232 

0.  0  0  03 

12.6000 

0.C250 

0.0006 

16.0000 

0.0270 

0.000c 

20.0000 

0. C298 

C.  0  0C8 

<x 


-4. OOCO 

O.COCC 

4.0000 

5.20CC 

8.00CG 

10.0000 

12.60CC 

16.0000 

20.0000 


c 

op. 

0, 

.0045 

0, 

COOO 

0. 

.0045 

0, 

C055 

0. 

C110 

0. 

0145 

0. 

C195 

0. 

0250 

0. 

C320 

~~t)  <* 


0.0011 
0.0  011 
C.00C9 
0.0  012 
C.  0018 
0.0  018 
0.0  018 
0.0017 
0.0018 
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TABLE   XV 
Correspondence  of  Symbols    (   nonlinear   Modal   ) 


X 
I 

q 

e 

Y1 
12 
X1 
X2 

tv 


X1 
X2 
X3 

xa 

X5 

X6 

X7 

X8 

X9 

X10 

X11 

X12 

X13 

X1U 

X15 

X16 

X17 


P 

cnsy 


rY*V 

C^y 

> 

?V 

CY(* 

qSd/IXx 

qSd/IZ2 

qSd/W 

qSd/IYy 

3.1416 

X18 
X19 
A1 

A2 

A3 

AU 

A5 

A6 

A7 

A8 

3A 

A  10 

CO 

C01 

C02 

C1 

C2 
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TABLE  XVI 

Correspondence    of  Symbols  (Cont.    Nonlinear   Hodel) 

188.4  C3  -~ =sj DA16 

150-0            C4  rN 

6.0  C5  ^-SmiSftZ      DA1a 


0.530544 C6 


TU 


r 


0.48  C7  ^K.— f-C-       DA26 


0.38375       C8 


t .. 


n> 


e*. 


3.07  C9 DA24 

5.0  C10  r\^^ 

0.05  033  C11  --- ***•--'-       DA3 

0.755  C12  rv     u 

8.0  C13  B«£ DAu 

17.6  C14  r-N         ^ 

6.0  C15  &&€*&'-      DA5 

0.078  C16  ^ 

15.0  C17  -— ^>-X.~ DA6 

4.0 C18  o,?0*    , 

0.913  C19  — orr—       DA7 

0.25  T1  pv 

0.839  K1  --^--P DA8 

6.08  K2  r\~       r    * 

1.0  KYP  —  ---£ DA9 


DAA 
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